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ABSTRACT 


The  effectiveness  of  the  use  of  front-end  bandstop  filters 
in  coherent  digital  receivers  versus  previously  derived  optimum 
jammer  waveforms  is  analyzed.   The  receivers  studied  are  the 
binary  frequency  shift  key  (FSK)  and  binary  phase  shift  key 
(PSK)  coherent  (correlator)  receivers.   The  filters  analyzed 
for  use  against  an  optimum  jammer  are  a  single  bandstop  region 
ideal  filter,  and  a  second  order  single  bandstop  region  real 
filter.   General  expressions  in  frequency  domain  form,  are 
derived  for  the  probability  of  error  of  coherent  receivers 
with  a  front-end  filter,  then  applied  specifically  to  the 
performance  of  receivers  operating  in  the  presence  of  noise, 
jamming,  and  PSK  or  FSK  modulation.   Finally  these  results  are 
presented  in  graphical  form  and,  additionally,  analyzed  and 
interpreted. 
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I.   INTRODUCTION 

This  thesis  investigates  the  use  of  filtering  techniques 
as  an  electronic  counter  countermeasure  (ECCM)  to  overcome  the 
effects  of  the  optimum  jammer  operating  against  coherent  phase 
shift  key  (PSK)  and  frequency  shift  key  (FSK)  receivers.   These 
receivers  are  well-known  solutions  to  the  problem  of  optimum 
detection  of  known  signals  in  the  presence  of  additive  white 
Gaussian  noise  (WGN)  [Ref.  1] .   The  structure  of  such  receivers 
is  shown  in  Figure  A.l.   These  two  receivers  represent  the 
optimum  signal  processing  algorithms  for  the  discrimination 
of  the  digital  signals  s, (t)  and  sfi(t)  representing  digital 
"one"  and  "zero"  logical  states  respectively.   To  reduce  the 
effectiveness  of  these  optimum  receivers,  Ref.  2  and  Ref.  3 
developed  optimum  jammers  to  be  used  against  the  PSK  and  FSK 
coherent  receivers.   From  the  resultant  jammer  waveforms  derived 
in  these  references,  this  thesis  investigates  the  use  of  front- 
end  filters  inserted  in  these  receivers  as  shown  in  Figure  A. 2 
as  a  method  for  reducing  jammer  effectiveness. 

From  Ref.  2,  the  optimum  power  constrained  jammer  for 
either  the  coherent  PSK  or  FSK  receiver  is  a  signal  propor- 
tional to  the  correlator  signal,  shown  in  Figure  A. 2  as  s^(t) . 
This  correlator  signal  (see  Ref.  1)  is  proportional  to  the 
difference  of  s, (t)  and  sn(t) ,  the  digital  "one"  and  "zero" 
signals.   Though  the  correlator  signal  can  take  different  forms 
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in  these  receivers  depending  on  the  modulation  used,  the  theory 
used  to  develop  these  receivers  shows  that  a  jammer  waveform, 
proportional  to  the  difference  of  s,(t)  and  sQ(t),  is  optimal. 

The  jammer,  therefore,  takes  the  form  of  a  continuous  wave 
(cw)  or  tone  centered  at  the  frequencies  of  s, (t)  and  sfl (t) . 
In  the  case  of  the  PSK  coherent  receiver,  the  jammer  is  a  single 
tone  centered  at  the  signaling  frequency.   For  the  FSK  case, 
the  jammer  waveform  consists  of  the  difference  of  two  tones 
at  the  frequencies  of  s,(t)  and  sQ(t).   The  use  of  tones  or 
cw  jammers  in  angle  modulated  analog  or  digital  receivers  is 
discussed  in  Ref.  4. 

Given  these  forms  of  jammer  signals,  an  obvious  choice  for 
the  front-end  filter  is  a  band  reject  or  bandstop  filter.   The 
filters  chosen  for  analysis  in  this  thesis  are  the  ideal  band- 
stop  filter  and  a  single  zero  second  order  bandstop  filter. 
Though  the  ideal  filter  is  not  realizable,  the  analysis  involv- 
ing its  use  will  provide  both  insight  into  the  problem  and 
approximate  results  before  investigating  the  use  of  the  more 
complex  second  order  filter.   For  the  PSK  coherent  receiver, 
the  bandstop  regions  of  both  the  ideal  and  second  order  filters 
are  centered  at  the  frequency  of  s, (t)  which  is  also  equal  to 
the  frequency  of  sn  (t)  .   Since  in  this  thesis  a  filter  with 
only  a  single  bandstop  region  will  be  analyzed,  the  bandstop 
region  of  the  ideal  and  second  order  filters  used  in  the  FSK 
coherent  receiver  will  be  centered  midway  between  the  frequen- 
cies of  s, (t)  and  sn(t) . 
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The  analysis  of  the  two  forms  will  begin  with  the  statis- 
tics of  G,  the  output  of  the  coherent  receiver.   From  Ref.  2, 
G  is  a  conditional  Gaussian  random  variable  with  conditional 
statistics  depending  on  whether  s, (t)  or  s~ (t)  was  transmitted. 
The  analysis  will  then  focus  on  determining  the  statistics  of 
the  random  variable  G1,  the  output  of  the  receiver  with  a  front- 
end  filter  as  shown  in  Figure  A. 2.   The  mean  and  variance, 
conditioned  on  whether  s, (t)  or  sn(t)  was  transmitted  will 
then  be  determined  using  frequency  domain  analysis  techniques. 
Finally,  the  general  analysis  will  conclude  with  the  determina- 
tion of  an  expression  for  the  probability  of  error  of  the  coherent 
digital  receiver  using  a  front-end  filter  to  counter  the  opti- 
mum jammer.   Next,  the  specific  frequency  domain  forms  of  the 
PSK  and  FSK  signals  and  jammers  will  be  developed  and  used  to 
calculate  the  probability  of  error  performance  for  these  re- 
ceivers with  and  without  front-end  filters. 

To  determine  the  Fourier  transforms  of  the  jammers,  this 

thesis  will  first  assume  that  the  jammer  waveform  is  a  time 

truncated  function.   Since  in  References  1,  2,  and  3,  the 

correlator  signal  s,(t)  is  defined  over  the  interval  (0,T), 

d 

and  the  optimum  jammer  (from  Ref.  3)  is  proportional  to  s,(t) , 
the  optimum  jammer  will  be  treated  as  defined  in  the  interval 
(0,T)  also.   Since  from  a  practical  standpoint,  the  jammer  is 
a  cw  signal,  the  analysis  will  be  modified  in  the  sequel  to 
account  for  this  interpretation  of  the  jammer. 

In  Chapter  II,  a  development  of  the  mathematical  expressions 
in  the  form  of  frequency  domain  equations  for  the  statistics 
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of  G'  (the  output  of  the  coherent  receiver  of  Figure  A. 2) 
is  undertaken.   From  these  statistics,  the  expression  for 
the  receiver  probability  of  error  is  developed.   Chapter  III 
demonstrates  the  application  of  these  general  expressions  to 
the  binary  PSK  (BPSK)  coherent  receiver  problem  with  both 
ideal  and  second  order  front-end  filters.   Chapter  IV  presents 
similar  analysis  to  that  in  Chapter  III,  applied  now  to  the 
binary  FSK  (BFSK)  coherent  receiver.   Chapter  V  presents 
quantitative  and  graphical  analysis  of  the  four  receiver/ 
filter  combinations,  under  the  assumptions  of  both  the  cw  and 
time  truncated  jammer  forms. 
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II.   ANALYSIS  OF  DIGITAL  COHERENT  RECEIVERS 
WITH  FRONT  END  FILTERS 

In  order  to  investigate  the  effects  of  filtering  at  the 
front-end  of  digital  coherent  receivers,  it  is  necessary  to 
review  the  operation  and  performance  of  digital  coherent 
receivers  in  the  presence  of  noise  and  jamming  without  the 
use  of  filters  to  notch  out  the  jammer.   As  previously  pre- 
sented in  Ref.  2,  the  structure  of  the  receiver  to  be  analyzed 
is  shown  in  Figure  A.l.   In  the  absence  of  a  jamming  signal 
this  receiver  is  optimum  for  determining  whether  a  "one" 
signal  or  a  "zero"  signal  was  transmitted  in  a  given  interval 
(0,T),  with  minimum  probability  of  error  [Ref.  1]. 

Following  the  analysis  of  Ref.  3,  which  assumes  that  the 
jammer  is  a  deterministic  waveform  unknown  to  the  receiver, 
the  input  signal  r(t)  to  the  receiver  front-end  is  mathematically 
modeled  as 


r(t)  =  s.(t)  +  n(t)  +  n.(t)   0  <_   t  <_   T,   i  =  0,1   (2.1) 


where  the  s. (t)  (i  =  0,1)  are  used  to  transmit  the  digital  "one" 
and  "zero"  data,  n(t)  is  a  sample  function  of  a  white  Gaussian 
noise  process  with  power  spectral  density  level  of  Nn/2  watts/ 
hertz,  and  as  stated  previously  n.(t)  is  the  jammer  waveform. 
From  the  analysis  of  Ref.  2,  the  output  of  the  receiver  is 
given  by 
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G  =   (r,sd)  +  |[| |sQ| |2  -  | |Sl| |2]  (2.2) 


where 


T  9 

(x,y)   =    /   x(t)y(t)dt   and    |x|T   =   (x,x] 


Since  G  is  a  (conditional)  Gaussian  random  variable,  the  sta- 
tistics are  completely  determined  when  the  conditional  mean 
and  variance  are  found.   The  conditional  mean  is  given  by 

E  [G  I  s  .  transmitted]   =   (s.,s,)  +  (n.,s,) 
1  1  id      j   d 

+  ^l|soM2-||s1||2]  (2.3) 

i  =  0,1 

and  the  conditional  variance  is  given  by 


N0        2 
VAR[G  |  s.  (t)  transmitted]   =   -y  ||sJ|  ,   i  =  0,1  .     (2.4 


A.   FILTERING  EFFECTS  ON  THE  MEAN  AND  VARIANCE 

The  effect  of  placing  a  filter  at  the  front-end  of  the 

receiver  is  now  analyzed  by  evaluating  the  conditional  mean 

and  variance  of  the  receiver  output.   Figure  A. 2  shows  the 

placement  of  the  filter  at  the  receiver  front-end.   It  can  be 
seen  that  G  becomes  G'  where 
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G-      =      (r',sd)    +  i[| [3o! |2  -  | |Sl| |2a 

(2.5) 

CO 

r'      =      r  *h      =  /      h(t-a)r(a)da 

—  oo 

and  h(t)  is  the  front-end  filter  impulse  response,   Since 


r'(t)   =   s!(t)  +  n'(t)  +  n!(t)   i  =  0,1      (2.6) 


it  is  simple  to  see  that 


E[G'|s.(t)  transmitted]   =   (s!  ,s..)  +  (n!,s,) 
l  i   d      j      d 


+  ji\   |b  | |2  - | |Sl| |2]      (2.7) 


i  =  0,1 


where 


s|(t)   =   s±(t)  *h(t) 


n'(t)   =   n.  (t)  *h(t) 


It  is  much  more  convenient  to  specify  filters  in  the  fre- 
quency domain  because  convolutions  in  the  time  domain  become 
multiplications  in  the  frequency  domain.   Thus  looking 
specifically  at  the  term  (s|,sd)  where 


T 

(s!,s,)   =    /   s!(t)s,(t)dt     i  =  0,1  (2.8) 

id      Q>         i    d 
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or 


(s^sd)   =    /   s£(t)sd  (t)dt    i  =  0,1         (2.9) 


where 


sdp(t)   =   sd(t)  •  p(t) 


with 


(1,       0  <_   t  <_   T 

P(t)   =   \ 

'  0 ,   elsewhere 


and  replacing   s,  (t)  by  its  inverse  Fourier  transform  equiva- 
lent (s,  (t)  <=>  S   (go)),  namely 


oo 

;,    (t)       =      -i-        /      S^    (uOe^dw  (2.10) 

dp  2tt      J         dp 


yields 


OO  00 

(s!,sj       =  /      s!(t)[^        /      S.    (co)e:|a3tda)]dt  (2.11) 

id  '         l  2tt         }         dp 


Interchanging  the  order  of  integration  yields 
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00  oo 

(s!,sd)  =2^   /   Sdp(aj)  [   /   s!(t)e:a)tdt]dw   i=0,l      (2.12) 


and  recognizing  the  fact  that 


/      s|(t)eja)tdt      =      S!  (-(d)  i   =    0,1 

— oo 


results    in 


1 
(s!,S,)       =      ±        j      S,    (a>)S!(-u>)du>         i   =   0,l  (2.13) 

i      a  ztt        '         dp  i 


Furthermore,  since 


sj(t)   =   si(t)  *  h(t)     i  =  0,1 


Sj(oj)  is  given  by 


S[(u)       =      Si(a))H(a))     i  =  0,1 


and   therefore 


(s!,s,)       =      -L        |       s       (o,)s.  (-w)H(-a))du)  i    =    0,1  (2.14) 

id  ziT         J  dp  1 
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Replacing  sd(t)  with  s,  (t)  as  in  Equation  2.9,  the  term 

(n!  ,  s  ,)  becomes 
j   d 


(n!,s,)   =    /   n!(t)s   (t)dt  (2.15) 

•J  _oo        J  " 


and  using  Equation  2.10  in  Equation  2.15  yields 


oo  oo 

(n!,sj   =    /   n!(t)[^-   /   S,  (u>)  e:ajtdu)]  dt  (2.16) 

3     d  -co        ^        27T   -J       dP 


As  done  in  Equations  2.11-2.13,  interchanging  the  order  of 
integration  and  using 


N!(-oj)   =     /   n!(t)eju)tdt  (2.17) 


with 


N!(oj)   =   N.  (u))H(o): 

D  D 


finally  yields 


(n!,s,)   =  ~   /   S,  (u)N.(-a))H(-a))du  (2.18) 

D   d      2tt  _m;    dp    : 
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The  conditional  mean  of  g',  given  that  s. (t)  is  transmitted, 
is 


E[Gi|s.(t)  transmitted]  =      ^-        \      S,  (a>)s  •  (-a)  H(-u>)du> 

1     l  2  it         *  dp  l 


+    o^T         /      S,(u))N.  (-<*>)  H(-(o)du> 


+   j[||soM2-||s1||2]  (2.19) 


i    =    0,1 


The  variance  of  g'  given  that  s. (t)  is  transmitted  can  be 
shown  to  take  on  the  form 


VAR[G'|si  transmitted]   =   E [ (n ' ( t) , sd( t) ) 2]  (2.20) 


The  right  hand  side  of  this  equation  becomes 


E[(n' (t) ,sd(t))2] 


=   E[   /   n'(t)sd  (t)dt    /   n'(T)sd  (x)dx]  (2.21) 


Equation  2.21  can  now  be  written  as 


E[(n' (t) ,sd(t))2] 


00  00 


/  /      E[n' (t)n' (t) ]sd    (t)sd    ( t ) dtdx  (2.22) 


—  OO  —  CO 
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Recognizing  that  the  expected  value  operation  in  the  integral 
is  the  autocorrelation  function  of  the  filtered  white  Gaussian 
noise,  R  .(t-x),  Equation  2.21  can  be  written  as 


E[(n\sd)2] 

OO  00 

/  /      Rn, (t-T)sd    (t)sd    (T)dtdx  (2.23) 


Again,  it  is  much  more  convenient  to  continue  the  analysis 
in  the  frequency  domain.   Using  the  inverse  Fourier  transform 
for  the  frequency  domain  representation  of  R  , (t-f)  and  sub- 


stituting into  Equation  2.23  yields 


E[(n\sd)2] 


00         00 


/     /  [^    /   Sn,(aJ)e^(t-T)da)]s   (t)s   (T)dtdx 


—oo     — oo 


dp     dp 

(2.24) 


where  s,  (t)  is  defined  by  Equation  2.9.   Interchanging  the 
order  of  integration,  Equation  2.24  becomes 


E[(n\sd)2] 

oo  oo  .  oo 

"      2¥        /     Sn'(w)dw         /      sdpltlell"tdt         /      sip^)s~^T^ 

— oo  — oo  "  —oo  c 

(2.25) 
It  is  easily  recognized  that 
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/   Wt)e3"tdt  -  Wa)   =  sdP<w) 


and 


/   s,^(T)e'Da)TdT   =   S,  (go) 
/         dp  dp 


Therefore  in  the  frequency  domain,  VAR(G'|s.  transmitted) 
becomes 


00 

VAR[G'|s-(t)  transmitted]   =   ~    /   S  ,(go)|s,  (go)  I  2dG0    (2.26) 
1  2tt   J         n     '  dp    ' 


Finally  since  n'(t)  =  n(t)  *  h(t),  it  is  well  known  that 


Sn,(oo)   =   Sn(o3)  |H(u)  |2 


Since  n(t)  is  a  sample  function  of  a  random  process  that  was 
assumed  to  have  power  spectral  density  S  (go)  =  N  /2,  for  all  go, 


Var[G'|s.(t)  transmitted]   =   ^   /  |h(go)|2|s,  (co)  |  2doo     (2.27) 
1  l  2tt  -.'  '     '  '  dp    ' 


2         2 
where  the  fact  that  |h(go)  |  "|S  ,  (go)  |   is  an  even  symmetric 

function  has  been  used  to  reduce  the  limits  of  integration 

to  half  the  real  line. 

22 


In  summary,  G'  is  a  (conditional)  Gaussian  random  variable 
with  conditional  mean  and  variance  given  by 


E[G'|si(t)  transmitted] 


-   2?   •/   Sd  (a))Si(-a))H(-co)doo 


+    2T         /       S       (a,)N    (-o))H(-ui)do)  +|[|  |s    |  |2  -  |  |s     |  |2] 

— oo         "       -J 

i  =  0,1  (2.28) 


and 


VAR[G'|s.(t)  transmitted] 


N     °° 
0   r   i 2,„   ,  %  i  2 


=   27  rt/   |H(o>)  r|Sd  (a))  rdu     i  =  0,1  (2.29) 


0 


B.   PROBABILITY  OF  ERROR  PERFORMANCE 

The  analysis  thus  far  has  been  concerned  with  determining 
the  statistics  of  the  random  variable  G1,  the  output  of  the 
coherent  receiver.   The  next  step  is  to  determine  the  proba- 
bility of  error  from  the  decision  process,  where  the  amplitude 
of  G1  is  compared  with  a  threshold  y.   Referring  to  Figure  A. 2, 
if  the  value  of  G1  is  greater  than  y,    the  receiver  decides  that 
a  "one"  or  s, (t)  was  transmitted,  and  conversely  if  the  value 
of  G1  is  less  than  y    the  receiver  decides  that  a  "zero"  or 
Sq (t)  was  transmitted.   The  quantity,  y,    is  determined  in  the 
derivation  of  the  optimum  receiver  for  detecting  known  signals 
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in  additive  white  noise  [see  Ref.  1]  .   That  analysis  determines 
the  threshold  of  the  correlator  (coherent)  receiver  as 


N  T 

Y   =^-£nri+j/  [y* (t)  -  y* (t) ]dt  (2.30) 

However  the  receiver  of  Figure  A. 2  treats  the  second  term  of 
Equation  2.30  as  a  bias,  thus  the  threshold  becomes 


N0 
Y   =  ~y   in   n  (2.31) 


where  N./2  is  the  power  spectral  density  level  of  the  white 
Gaussian  noise  process.   Here,  r\    is  determined  by  the  specific 
decision  rule  applied.   In  this  case, the  receiver  strategy 
is  to  detect  s, (t)  or  sn(t)  with  minimum  probability  of  error, 
so  from  Ref .  1 


P{s, (t)  was  transmitted} 
P{sn (t)  was  transmitted} 


All  further  analysis  in  this  thesis  will  assume  that  the 
probability  that  s, (t)  was  transmitted  is  equal  to  that  of 
sn (t)  being  transmitted,  so  the  decision  threshold  is 

N0 
Y   =  -V.   £n  1   =   0  (2.33) 

Now,  the  receiver  can  make  two  types  of  errors.   Specif- 
ically, deciding  that  sQ(t)  was  transmitted  when  in  fact  s-^t) 
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was  transmitted  (defined  as  error  1) ,  and  deciding  that  s, (t) 
was  transmitted  when  sQ(t)  was  actually  transmitted  (defined 
as  error  2).   The  receiver  decides  that  s~ (t)  was  transmitted 
when  G'  is  less  than  the  threshold  y,  and  that  s, (t)  was  trans- 
mitted when  G1  is  greater  than  y.      So,  the  probability  of  error 
1  occurring  is  the  probability  that  G1  is  less  than  y   when  s, (t) 
is  transmitted.   Since  G'  is  a  Gaussian  random  variable  with 
conditional  probability  desnity  functions  dependent  on  whether 
s, (t)  or  sn(t)  was  transmitted 

(g'-rn^2 

Y    i         2v 

P      .   =     /     X    e      g    dg1  (2.34) 

error  1       >        ,-* r 

g 


where  m  =  E[G'|s,(t)  was  transmitted]  and  v   =  VAR[G'|s1(t) 
was  transmitted].   By  letting  x  =  g'-m,//v  ,  Equation  2.34 
becomes  the  integral  of  the  standard  Gaussian  probability 
density  function  also  known  as  the  error  function  (erf(x)) 
of  the  form 


-m1//vg      ^  2 

P       ,   =     f         -±-   e  X  /2  dx  (2.35) 

error  1       ;         ,-*— 

-oo  /2tt 


The  probability  of  error  2  occurring  is  the  probability  that 
G'  is  greater  than  y   when  s  (t)  is  transmitted.   Therefore 

2 
i    -(g'-mn)  /2v 

P       ,   =    f  —L—   e       °      g  dg'  (2.36) 

error  2       *         ,~ 

y  /ZfTV 

g 
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where 


nu   =   E[G'|s.(t)  transmitted]. 


Using  a  similar  change  of  variables  as  in  Equation  2.35 
yields 


oo  2 

P       0   =         /   -=-  e"X  /2  dx  (2.37) 

error  2         ,  , '   pr— 

0    g 


where  y  is  zero  from  Equation  2.33.   Finally,  the  average 
probability  of  error  becomes 

P    =   P  P       ,    +    d-p)  P       o  (2.38) 

e     v      error  1      v        error  2 

where  p  =  P{s, (t)  was  transmitted}  and  1-p  =  P{s0(t)  was 
transmitted}.   Since  p  =  1/2  and  substituting  Equations  2.35 
and  2.37  into  2.38,  the  average  probability  of  error  is 

-m,  //v~ 


1    g   '     -z 

e  2 


P.     =     ki      I      *        y^e-'/2dx 


-oo  /2tt 


-mn//v  /2tt 

0         g 


2 
J      -i-  e"X   /2    dx]  (2.39) 
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III.   PHASE  SHIFT  KEY  COHERENT  RECEIVER 
WITH  A  FRONT-END  FILTER 

A.   PHASE  SHIFT  KEY  WITH  AN  IDEAL  FRONT-END  FILTER 

Using  the  equations  for  the  mean  and  variance  of  the  con- 
ditional Gaussian  random  variable  G*  (Equations  2.28  and  2.29), 
along  with  the  probability  of  error  equation  (Equation  2.39) , 
it  is  possible  to  investigate  the  effect  of  a  front-end  filter 
on  the  probability  of  error  performance  of  the  coherent  receiver 
The  receiver  analyzed  in  this  case  is  the  optimum  receiver  for 
binary  phase  shift  key  modulated  signals  with  an  ideal  band- 
stop  filter  in  the  front-end.   The  filter  will  be  characterized 
as  having  a  gain  of  unity  over  all  frequencies,  with  the  excep- 
tion of  the  regions  from  cu  =  ±(jon  -B/2  to  oj  =  ±oo0  +B/2  where 
the  gain  is  zero.   The  bandwidth  of  the  filter  is  B  radians/ 
second  and  the  center  frequency  is  con  radians/second.   Although 
this  filter  is  not  realizable,  it  serves  both  as  a  good  approxi- 
mation to  a  higher  order  bandstop  filter,  and  as  a  simple 
example  so  as  to  gain  insight  to  the  receiver  performance. 
Using  Figure  A. 2  as  a  reference,  the  generic  coherent 
digital  receiver  can  be  made  into  a  BPSK  receiver  by  setting 


s,(t)   =   2A  cos  coAt     0  <  t  <  T  (3.1) 

d  0        —   — 


and  modeling  the  input  PSK  signals  as 
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si(t)       =      A  cos(ooQt   +    9  .  (t)  )       0    £  t    <    T,       1   =    0,1  0.2: 


where 


fir      for      i   =   0 
0      for      i   =    1 


Thus    s.(t)     (i   =    0,1)    can  be    expressed   as 


s0(t)       =      -A   cos   co^t 


s,  (t)       =      A   cos   u)fit 


The  jammer  waveform  is  set  to 


0    <    t    <    T  (3.3) 


n.(t)       =      /P      '   /2/T   cos   u>At         0    <    t    <    T  (3.4) 

j  nj        '  0  —        — 


since  in  Ref.  2  it  was  shown  to  be  optimum  against  BPSK  trans- 
mission.  Here  P    is  the  jammer  output  power.   Finally  the 
noise  input  into  the  receiver  of  Figure  A. 2  is  a  white 
Gaussian  random  process  with  power  spectral  density  level  of 
N^/2  watts/hertz. 

1 .   Calculation  of  the  Variance 

From  Equation  2.27,  the  variance  of  G'  conditioned  on 
s.(t)  being  transmitted  is  given  by 
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VAR[G'|s.(t)    transmitted] 


N«        °° 


0       f       i  „  ,    %  ,  2  ,  „       ,    ,  ,  2 


=      «-     /       |h(u>)  |     |S-    (u)  |    du        i   =   0,1  (3.5) 

Z7r    n  aP 


0 


Since  the  front-end  filter  is  ideal,  Equation  3.5  can  be 
modified  to  become 


VAR[Gf|s.(t)  transmitted] 

N  -  u)0+B/2 

=      57 [    /     |S       (oo)  |2dco    -  /  |     S       U)|2du3  (3.6) 

0  P  oo0-B/2  ap 

i  =  0,1 

2 

since  |H(oj)  |   is  unity  except  over  the  bandstop  region. 

i  2 
In  order  to  calculate   S,  (a))   ,  since  s.  (t)  =  s,(t)p(t) 

1  dp    '  dp       d  '  r      > 

where 


!1   0  <_   t  <  T 
0   elsewhere 


then 


sd  (t)   =   2A  p(t)  cos  pQt  (3.7) 


The  Fourier  transform  identity   f(t)-g(t)  <=>  (l/2ir)  xf(co)  *  G(co) 

will  be  used  in  order  to  calculate  S,  (co)  .   Since  P(co)  is 

dp 
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PU)       =      T   Sa(ooT/2)e   ^T/2  (3>g) 


and 


Sd(co)       =      2Att  [6  (oo-oo0)    +    6(o)+co0)J 


S,  (oo)  is  the  convolution  of  these  Fourier  transforms,  result- 
dp 

ing  in 


S   (w)   =  AT[Sa(a)-coQ)^e 


-j  (oo-oon)T/2 


T  -j  (a)+oo0)T/2 


+  Sa  (w+a)0)  je  ]  (3.9) 


2 
From  this,  |  S  ,  (co)  |   is  obtained  by  multiplying  S,  (u))  by 

S,  (-03)  which  results  in 
dp 

|Sd  (oo)  |  2   =   A2T2  •  [Sa^oo-WQ)  ^ +Sa2(o)+oa0)|]  (3.10) 


where  the  product  of  Sa  (oa+ca.)  T/2  and  Sa  (oo-oo.)  T/2  has  been 
ignored  because  it  is  essentially  zero  for  practical  values 
of  u)n  and  T. 

Since  the  variance  consists  of  an  integral  from  oo  =  0 
to  infinity,  for  large  oon  the  contribution  of  the  Sa  (oj+w^)  T/2 
term  is  negligible.   Therefore,  substituting  into  Equation  3.6 
yields 
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Var [G ' | s . (t)    transmitted] 


2    2   Nn        °°       ?  t  °  5  T 

=    A   T      ^-[    f      Sa(a)-a)n)^du)  /  Saloo-con)^u)]  (3.11) 

27T       f\J  0       2  T3/0  0       2 

0  0)q-B/2 

i    =    0,1 

In  this  form  Equation  2.11  is  unmanageable  due  to  its  complexity 
and  the  number  of  variables  present.   To  simplify,  the  follow- 
ing substitutions  will  be  made: 

1.  x  =  (o)-oo0)T/2/   2/Tdx  =  doo 

2.  for   03    =    0 ,    x   =   -conT/2   which    for    large   values   of   oonT 
becomes    approximately    -co 

3  .       for  oo  =  co  ,    x   =  co 


4.   for  oo  =  oo   ±B/2,  x  =  ±BT/4 


Equation  3.11  now  becomes 


VAR[G'|s.(t)  transmitted] 


A  TNn     co  BT/4    2 

=   -[      j      Salx)dx  -       /      Sa(x)dx]  (3.12) 

77   -co  -BT/4 

i  =  0,1 

2 
Finally  the  Sa(x)  is  a  symmetric  function  about  x  =  0  and 

°°   2 
|  Sa  (x)dx  =  tt/2  from  math  tables.   Thus  the  variance  of  G1 

0 


31 


can  be  written  as 


VAR[G'|s.(t)  transmitted] 


BT/4 
=   N  AZT  -[y  -   /     Sa^(x)dx]    i  =  0,l  (3.13) 

U       TT  Z.  q 


2 .   Calculation  of  the  Conditional  Mean 

The  conditional  mean  of  G'  is  given  by  Equation  2.7 
and  shown  here  for  convenience, 


E[G*|s.(t)  transmitted] 


=   (s!,sd)  +  (n!,sd)  +  |[|  |s0|  |2  -  |  |Sl|  |2]    i  =  0,1 

and  consists  of  three  terms.   The  frequency  domain  forms  of 

(s!,s,)  and  (n!,s,)  are  given  by  Equations  2.14  and  2.18 

12         2 
respectively.   The  term  j[\   |  s~  |    -  |  | s,  |   ]   (recall  that 

2  5 

| sn |    =   /   sn(t)dt),  represents  the  difference  between  the 
u      0 

energy  contained  in  sn(t)  and  s, (t) .   Since  our  signals  are 
equal  energy  signals,  the  term  in  question  becomes  zero.   Thus 
the  conditional  mean  becomes 


E[G'|s.(t)  transmitted]   =   (s!,sd)  +  (n!,sd)   i  =  0,1     (3.14) 


Looking  first  at  the  term  (n!,s,)  and  recalling  that 

J 

H(w)  is  an  ideal  filter,  Equation  2.18  can  be  modified  to 
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(ni'Sd}       =      2¥[       /      S       <u>)N    <-u,)da> 


oa0+B/2  -oo0+B/2 

/  S,     (uj)N.  (-co)doo  -  /  S,     (go)N.  (-a>)du>: 

^0-B/2;  d?  J  -u)0-B/2J  dP  ^ 

(3.15) 


In  order  to  evaluate  Equation  3.15,  the  product  N.(-u))S,  (u>) 
must  be  calculated  using  n.(t)  as  given  by  Equation  3.4. 
However,  since  the  model  of  the  coherent  receiver  in  Figure 
A. 2  contains  a  finite  time  integrator,  n. (t)  is  effectively- 
truncated  to  the  interval  0  <  t  <  T,  or  equivalently ,  the 

new  term  n.  (t)  replaces  n.(t)  where  n.  (t)  =  n.(t)p(t)  and 
3P  3  DP        3         * 

p(t)  is  as  previously  defined.   Thus,  N.  (co)  becomes 

J  kr 


N.     (oj)       =       /P         /2/T    i- 
DP  n:  2 

-j  (go-go    )T/2  -j  (uW-co    )T/2 

•     [Satoo-w^^e  +    Sa(co+co0)T/2e  ]  (3.16) 


Multiplying  N.  (-co)  and  S,  (co)  yields  four  terms:   two  Sa  terms 

centered  at  oon  and  -oon  and  two  terms  that  involve  the  product 

of  Sa(co+oo  )T/2  and  Sa  (co-oj  _)  T/2  .   It  can  easily  be  seen  that 

for  large  ajn  this  product  is  essentially  zero.   Substituting 

into  Equation  3.15  it  is  also  easy  to  see  that  the  contribution 

of  the  Sa  (oo+oo^)  T/2  term  in  the  integral  with  limits  of  +co  +B/2 

is  negligible,  so  that  (n!,s,)  becomes 

J 
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1    AT2 


(nrsd'    -   2¥^/pnj  /2/T 


00 


o)Q+T/2 


2,  x    T  ,  (  „    2 ,  XT 


/£  X  (  £  X 

Sa    ((jo-oj-)   -jdoo    -  J  Sa    (co+(joq)   ydoo) 

co0-B/2 

"W0+B/2 
+  (   /  Sa2(o)+u0)  |-        /        Sa2(o3+oo0)  |dco)  ]       (3.17) 
-°°  -con-B/2 


2  2 

Recognizing  that  the  integrals  of  Sa  ((D+u)n)  T/2  and  Sa(d)-u-)  T/2 

are  equal,  performing  the  same  change  of  variables  as  in 

Equation  3.11,  (n!,s_J  becomes 


AT 


(n!, s,)   =   ~  /P  .  /2/T 
j   d       tt    nj 


oo  BT/4 

[   /  Sa2(x)dx  -       /     Sa2(x)dx]  (3.18) 

-oo  -BT/4 


Finally,  recognizing  symmetry  about  x  =  0  and  using 

0 


/  Sa  (x)  dx  =  tt/2,  we  obtain 


9        BT/4 

(n! ,s,)  =  AT  /P~7  /2/T  -[^  -   /     Sa^(x)dx]  (3.19) 

j   d  n:       tt  2    Q; 


Looking  now  at  the  term  (s!,s,),  recalling  again  that 
H(co)  is  an  ideal  filter,  (s!,s,)  can  be  written  as 
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■i     °°  0  ' 

(si/Sd)   =   2~¥[   /  sdp(w)Si(-a))dcJo  /       Sd  (aj)Si(a))doj 


-oa0-B/2 


J         S,  (u))S.  (-o))du)]  (3.20) 

-o)0-B/2J  dP 

i  =  0,1 


The  signal  models  for  s_(t)  and  s,  (t)  given  by  Equation  3.3 
can  be  Fourier  transformed  to  yield 


a  ,  ,         ATra  ,      ,t  -3<"-"0)T/2 

Sq  (u>)   =   -  ~2~  [Sa(co-cog)  ^e 

-j  (o)+o>  )T/2 

+  Sa(oo+oo0)  i-e       u     ]  (3.21) 


and 


AT  -j(o)-o)  )T/2 

S,  (oo)   =   •y[Sa(ai-(i)Je 

-j  (a)+co  )T/2 

+  Sa(oo+oo  )  i-e  (3.22) 


To  calculate  the  specific  case  of  (s!,s,),  first  form 
the  product  of  S.  (u))Sn(-u)).  This  yields  four  terms  as  in 
(n!,s,),  where  two  are  S  a  terms  centered  at  ioj-  and  two  are 
products  of  Sa  (oj+oon)  T/2  and  Sa  (w-co0  )  T/2  .   The  latter  products 
can  be  neglected  for  large  con .   Also,  with  respect  to  Equation 
3.20,  the  contribution  of  the  Sa (oo+con)  T/2  term  to  the  integral 
whose  limits  are  +ajn+B/2  is  negligible.   Therefore,  substituting 
S,  ((jo)Sn(-co)  into  Equation  3.20  yields 
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2    2  <=° 

(s0'Sd}       =      "A2T      2¥^  (      /    Sa2(oo-co0)T/2    doo 

—  oo 

wo+B/2 

/  Sa    (u>-U)0)-=do))    +    (       /  Sa  (oo+oa    )jdca 

00--B/2  -°° 

-a)0-B/2 

/  Sa2(co+oon)   Sdcu)  ]  (3.23) 

-w0-B/2  ° 


It  is  easily  seen  that  the  integrals  involving  Sa  (oj+u)n)  T/2 


2 
l( 

are  equal  to  the  integrals  involving  the  SdT(o3-U)0)  T/2  term. 

Making  use  of  this  observation  and  the  same  change  of  variables 

used  in  Equation  3.9,  it  is  easily  seen  that 


fo.  c,  )         -        -A2T  I  °°   9  BT/4   2 

^S0' Sd;  ir[   /  Salx)dx  -       /     Sa(x)dx]        (3.24) 

-BT/4 


2 

and  by  using  the  symmetry  of  the  Sa (x)  function,  (s',s,) 

becomes 


,  7  BT/4 

(s'sj   =   -AZT  -[^-  -   /     Sa^(x)dx]  (3.25) 

0    Cl  7T  l  ' 


By  similar  analysis  it  can  be  shown  that 


7  BT/4 

(s_\s,)       =      AZT   -[^   -       /  SaMx)dx]  (3.26) 

Id  it    2.  J 
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In  general  terms,  (s!,sd)  (i  =  0,1)  can  be  written  as 


BT  /4 
(s'sJ   =   -l(l+1)A2T  hi  -      j  Sa2(x)dx]   1  =  0,1      (3.27) 

0 


so  that  the  conditional  mean  of  G*  can  be  expressed  as 


E[G'|s.(t)  transmitted] 


RT  /4 

[AT/P-   /27T   +-l(l+1)A2T] [-{i  ~      I  Sa2(x)dx)]  (3.28) 

nj  7T    2         Q^ 

i   =    0,1 


3 .   Probability  of  Error  Calculation 

With  the  conditional  mean  and  variance  of  G1  given  by 
Equations  3.28  and  3.13,  it  is  now  possible  to  compute  the 
probability  of  error  performance  of  the  BPSK  coherent  receiver 
with  an  ideal  front-end  filter.   Assuming  the  probability  of 
the  transmitter  sending  a  "one"  or  a  "zero"  signal  is  the  same 
(P(s,(t))  =  P(sQ(t))  =  1/2)),  Equation  2.39  shows  that  the 
average  probability  of  error  P   is 


oo         2  -m,//v      2 

Pe  =  j[  j      —  e'x  /2dx  +    /        e"X  /2dx]     (3.29) 

-m./Zv-    /2tF  -oo 

0    g 


where  the  threshold  setting  given  by  Equation  2.33  has 
been  utilized. 
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In  order  to  compute  P  ,  all  that  is  needed  is  the 
specific  conditional  mean  divided  by  the  square  root  of  the 
conditional  variance  and  using  this  ratio  in  the  limit  of 
the  appropriate  error  function  integral.   First,  for  a  "zero" 
transmitted, 


9         0  BT/4 

[AT/Pn./27T-A^T]  [=■&-     I  Sa^(x)dx)] 

-mQ  nD  7T    2         > 

U (3.30) 


[  —  r- 

'0  7TL2 


-  /AZTNn  £[J-  /     Sa^(x)dx) 


0 


2 
Recognizing  that  A  T/2  =  E,  the  average  energy  per  bit  and 

substituting  this  in  Equation  3.30  one  obtains 


2  ,,      BT/4 
_m       [/4EP    -  2E] [~(|  -   /     Sa^(x)dx)] 


^  r  o  -n  BT/4 

0  ttv2 


(3.31) 


-  /2ENn  ±&  -      j  Sa     (x)dx) 


Finally 


-mn   L  BT/4 

— ^  /2E/Nn(l-/Pn./E)A/|-(y  -   /     Sa^(x)dx  (3.32) 

/v  J    V71  z    o 

g 


and  similarly 


m,  L  BT/4 

-     v2E/NQ(l    +   /Pn./E)W^-(j  -      /  Sa^  (x)dx)  (3.33) 


v  ~  J  v        -         o 

g 
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By  using  Equations  3.32  and  3.33  as  the  limits  of  integration 
in  Equation  3.29,  the  probability  of  error  performance  can  now 
be  determined  versus  E/N0 ,  the  signal  to  noise  ratio;  P  ./E, 
the  jammer  to  signal  ratio;  and  BT,  the  relationship  of  the 
filter  bandwidth  to  the  inverse  of  the  signal  bandwidth. 

Several  items  of  note  can  be  observed  in  the  calcula- 
tion of  P  .   First,  as  BT  approaches  zero,  the  term 

e   BT/4   2 
2/tt  ( tt/2  -   /     Sa  (x)dx)  becomes  unity  and  Equation  3.29 

0 
becomes  the  receiver  probability  of  error  under  optimum  jamming 

[See  Ref.  2]  for  BPSK  modulation.   Furthermore,  when  both  BT 
and  P  ./E,  the  jammer  to  signal  ratio,  go  to  zero,  Equation 
3.29  becomes  the  well-known  result  for  performance  of  a 
coherent  BPSK  receiver  operating  in  additive  white  Gaussian 
noise  [see  Ref.  1] . 

B.   COHERENT  PSK  USING  A  SECOND  ORDER  FRONT-END  FILTER 

Instead  of  choosing  a  filter  transfer  function  that 
approaches  the  ideal  case,  the  effect  of  a  second  order  front- 
end  filter  on  the  probability  of  error  performance  of  the 
coherent  BPSK  receiver  will  be  investigated.   The  general 
coherent  receiver  structure  becomes  the  BPSK  coherent  receiver 
by  using  Equations  3.1  and  3.3  in  place  of  s,  (t)  and  s.  (t) 
(i  =  0,1)  respectively.   The  filter  will  be  characterized  as 
having  a  single  notch  at  the  center  frequency  of  the  jammer. 
The  filter  transfer  function  H(s) ,  is  derived  from  a  first 
order,  low  pass  prototype  given  by 


H(s)   =  j±-  (3.34) 

1+s 
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To  transform  this  into  a  bandstop  filter  centered  at  con  ,  the 
frequency  transformation 


s(co   -  0)  ) 
s   =   -— jH —  (3.35) 

S   +  U  d)T 

U   Li 


will  be  substituted  in  Equation  3.34.   Here,  to   and  to   repre- 

U  Li 

sent  the  upper  and  lower  3  dB  frequencies  of  the  bandstop 
region.   Before  substitution,  it  can  be  seen  that  for  co 

L 

approximately  equal  to  oj  ,  or  equivalently  co  -co   small  com- 

2 
pared  to  co  to  ,  to  to   and  co  -  co   are  approximately  equal  to  co~ 

U  L    U  1j        U    Li  u 

and  B  respectively.   Note  that  B  is  the  half  power  bandwidth 
of  the  bandstop  region.   Thus,  Equation  3.3  5  becomes 


s   =   — 2 J  (3.36) 

S   +U3Q 

Substituting  finally  into  Equation  3.34  and  also  setting  s  =  jco 
in  order  to  arrive  at  the  frequency  domain  form,  the  filter 
transfer  function  is 

2    2 

CO   -OJ 

H(co)   =   -^ ^ (3.37) 

to.  -  co  +  jcoB 

This  filter  will  prove  somewhat  cumbersome  to  work  with  when 
put  into  Equations  2.14,  2.18  and  2.27.   Therefore  a  more 
convenient  form  of  this  filter  will  be  sought  later. 

Looking  at  Equations  2.14,  2.18,  and  2.27,  the  products 

2 

S,  (co)S.(-co),  S  .,  (co)N.(-co)  and  Is,  (co)  I   (calculated  in  the 
dp     l        dp     j  dp 
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derivation  of  Equations  3.10,  3.18,  3.23)  are  Foureir  trans- 
forms of  a  time  limited  signal  multiplied  by  a  cos(tont).   The 
time  domain  forms  of  these  products  are  of  the  general  form 


c(t)   =   a(t)cos(a)0t)     0  <_  t  <  T  (3.38) 


which  is  essentially  an  "AM"  modulated  signal.   This  is  not 
surprising  because  each  signal  in  the  coherent  BPSK  receiver 
model  is  of  constant  phase  over  a  bit  interval  (0  <  t  <  T) . 
Reference  5  shows  that  the  effect  of  filtering  a  modulated 
cosine  wave  can  be  determined  by  calculating  the  effect  of  the 
modulation  being  passed  through  an  equivalent  baseband  filter. 
The  baseband  equivalent  of  the  bandstop  filter  will  be  in  this 
case  a  high  pass  filter. 

The  equivalent  high  pass  filter  model  can  be  derived  by 
applying  the  transformation 


s   =   —  (3.39) 

s 


to  the  low  pass  prototype  of  Equation  3.34,  and  letting  0.    , 
the  upper  cutoff  frequency  be  equal  to  B/2.   The  resulting 
filter  in  the  frequency  domain  is 


H'(u)   =   B/23"  ju  (3-40» 


When  Equation  3.37  is  evaluated  under  the  assumptions  of 
large  con  and  small  B  compared  to  oon,  the  filter  of  Equation 
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3.40  can  be  seen  as  an  equivalent,  frequency  shifted  version 
of  the  filter  of  Equation  3.37.   Using  the  filter  of  Equation 
3.40  and  just  the  "modulation"  spectra  of  the  products,  the 
effect  of  a  second  order  bandstop  filter  in  the  front-end  of 
the  BPSK  receiver  can  be  evaluated  with  respect  to  the  proba- 
bility of  error  performance. 

1 .   Calculation  of  the  Conditional  Mean 

The  conditional  mean  of  G'  is  given  by  Equation  2.7 

where  (n!,s,)  and  (s!,sj  are  given  by  Equations  2.14  and 

2         2 
2.18,  with  the  added  fact  that  1/2 ( |  | sfi |    -  |  | s,  |  |  )  is  zero. 

There  is  now  enough  information  to  determine  the  conditional 

mean  of  G'  namely 


E[G'|s.  transmitted]   =   (n!,sd)  +  (s!,s,)   i  =  0,1 


Taking  each  term  individually,  (n!,s,)  will  be  calculated 
first.   Equation  2.18,  repeated  here  for  clarity, 


00 

(n!,s,)   =  ?f-         /   S,  (w)  N.  (-oj)H(-w)du) 
j   d      2tt    >         dp    2 


shows  that  the  product  S  ,  (oj)N.(-oo)  and  H(-co)  will  be  needed. 

dp     j 

The  product  S,  (u)N.(-a))  was  previously  found  (in  the  derivation 


of  Equation  3.17)  to  be 


2 
AT 


S,  (oo)N.(-u))   =   ^—  /P  ./2/TtSa^Go-u).)^ +  Sa  (oo+uj-)t]         (3.41) 
dp     j  2     nj  0  2         0  2 
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Before  multiplying  this  by  H(-oj)  and  substituting  into  Equa- 
tion 2.18,  it  can  be  seen  that  Equation  3.41  can  be  written 
in  the  form 


S,  (u>)N.  (-a)) 
dp    j 


2  ,    VT 


=   AT/P.  /2/T  TSa  (o))~   *   tt  [  6  (w-co0 )  +  6  (oo+u)Q)  ]  (3.42) 


which  in  the  time  domain  is  of  the  form  of  Equation  3.38. 
Thus  the  baseband  equivalents  of  S,  (co)N.(-co)  and  H(-w)  can  be 
substituted  into  Equation  2.18.   Therefore  (n!,s,)  becomes 


ou 

(n!,s-,)       =      ^-AT/P"- /27f         /   TSa2  (ooT/2)  H  ■  (-co)  dco  (3.43) 

i      d  2tt  ni  J  ' 


2 

where  TSa(coT/2)  is  the  modulation  spectrum  and  H'(-co)  is  the 

equivalent  high  pass  filter   Finally,  substituting  Equation 
3.40  into  Equation  3.43  yields 


(n!,s,)   =   -^AT/P  V2/T    /  TSa2  (coT/2)  •   .~^/0  dco       (3.44) 
j      d       2tt      nj      '  J  '  -;jco+B/2 

Now  taking  advantage  of  the  simplified  expressions  for  the 
filter,  the  "modulation"  term,  and  the  respective  time  domain 
forms,  Equation  3.44  can  be  evaluated  using  Parseval's  theorem. 
Parseval ' s  theorem  states  that  if  f. (t)  and  f?(t)  are  time 
functions  with  Fourier  transforms  F,  (co)  and  F~(co)  then 
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oo  oo 

/      f1(t)f2(t)dt      =     -L       /      F    (u)F2C-u)) 


dxo 


By   letting  -joi/ (- joj+B/2)    equal   F_  (-to)  ,    it    is    easy   to    show   the 
impulse   response   of   the  baseband   filter  h'(t)     is 


h'(t)       =      6(t)    -   |  e    Bt/2    u.(t)  (3.45) 


Before  using  Parseval * s  theorem,  the  inverse  Fourier  transform 

2 
of  T(Sa(to)T/2  must  be  found.   This  is  easily  determined  from 

tables  to  be 


|  +  1,   -T  <_   t  <_   0 
F"1[TSa^00T/2)  ]   =   j  (3.46) 

(  -  %  +  l,    0  <  t  <  T 


Putting  (n!,s,)  into  the  equivalent  time  domain  integral  form 
using  Parseval's  theorem,  yields 


T 
(n!,sd)   =   AT  /P~  /27T   /  (-£  +  !)  (6  ( t) -|e"Bt/2y  ( t)  )  dt    (3.47) 


due  to  the  fact  that  the  filter  is  causal  (only  defined  for 
t  >  0)  and  that  F_1  (TSa^u))  T/2)  is  defined  from  -T  <  t  <  T. 
Evaluating  Equation  3.47  consists  simply  of  the  task  of  evalu- 
ating the  sum  of  four  integrals  with  limits  from  zero  to  T. 
It  is  easily  demonstrated  that 
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-BT/2. 


(nj,sd)   =   AT  ^F~  /2/Tt^Cl  -e  Di/")]  (3.48) 


The  second  term  that  must  be  evaluated  is  (s!,s,), 

1   d 

(i  =  0,1)  in  order  to  be  able  to  determine  the  conditional  mean 
of  G'.   From  Equation  2.14  repeated  here  for  convenience, 


oo 

(s!,  s,)  =     ■£-        /   SJ  (a))S  .  (-oo)H(-oo)  doo   i  =  0,1 

1    a  2TT   -oo       dP       1 


It  is  apparent  that  the  products  of  S,  (oo)  and  S.(-co)  (i  =  0,1) 

dp         i 

have  been  determined  in  the  derivation  of  Equation  3.24  and 
can  easily  be  put  in  the  form  of  Equation  3.42.   Therefore, 
the  baseband  equivalents  of  the  filter  and  the  product  are 
substituted  into  Equation  2.14  which  yields 


^T      m»  2  ,,..„,,,%    -jh) 


(s0'Sd)   "   -%F   /  TSa^T/2)  ^3^. 

—  oo  ->  ' 


dw  (3.49) 


It  is  easily  seen  that  the  integral  in  this  equation  is  of 

the  same  form  as  Equation  3.44,  thus  Parseval's  theorem  can 

be  applied  directly  and  for  the  specific  case  of  (s',s. )  yields 


T 
(sn,s)   =   -A2T   /  (^|-+1)  (6(t)  -|e"BT/2y(t))dt  (3.50) 

U   d  0 


Here  the  limits  of  the  integral  are  0  and  T  because  the  filter 
impulse  response  h'(t)  is  causal.   Since  the  integral  has 
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already  been  evaluated,  we  obtain 


(s£,sd)   =   -A2T[^(1  -e"BT/2)]  (3.51) 


and  similarly 


(s{,sd)   =   A2T[^(1  -e  BT/2)]  (3.52) 


Finally,  the  conditional  mean  of  G'  under  the  assump- 
tion that  a  "one"  or  a  "zero"  signal  was  transmitted  is 


E [G '  s .  transmitted] 
1  1 

=       [AT    /P~   /ITT  +    -l"i  +  1    A2T]  [-^(1  -e"BT/2)  ]  (3.53) 

n  j  di 

i   =    0,1 


2 .   Calculation  of  the  Conditional  Variance 

In  order  to  calculate  the  conditional  variance  of  G' , 
it  will  be  more  convenient  to  use  Equation  2.27  due  to  the  fact 
that  Parseval's  theorem  can  be  applied  to  its  evaluation. 
The  conditional  variance  of  G'  is 


VAR[G'|s.(t)  transmitted] 

N      °° 
=   4F    /   |H  (oo)  |  2  |  S   (co)  |  2  dco    i  =  0,1  (3.54) 

-oo  P 

i       i  2  . 

The  term   Sn  (co)    has  been  calculated  and  is  given  by  Equation 
dp 

3.10,  and  as  shown  here,  it  is  of  a  form  similar  to  Equation 
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3.41,    namely 


Sdp(o))|2    =   A2T2[Sa^u)-w0)T/2+    Sa2  (oo+Wq  )  T/2]  (3.55) 


2 

Therefore  the  baseband  equivalent  of  |  S,  (to)  |  ,  that  is 


|S'  (oo)|2   =   2A2T2[Sa2(ooT/2)]  (3.56) 


will  be  used  along  with  the  baseband  equivalent  filter  H*(w). 

2 
|H'  (co)  I   is  calculated  by  multiplying  H'  (go)  by  H'  (-co)  ,  producing 

the  following  result 

?  2 

|H'(o))r   =   5 5"  (3.57) 

(B/2)Z  +  oo 

Thus  the  conditional  variance  of  G1  becomes 

N  2A2T 

VAR[G'|s.  transmitted]   =  -^— -. 

1  l  4tt 

o             2 
/  TSa^(o)T/2)[ ^ y]dw   i  =  0,1        (3.58) 

(B/2)^  +0) 

by  substituting  Equations  3.56  and  3.57  into  3.54. 

Use  of  Parseval's  theorem  in  order  to  compute  the 

conditional  variance  of  G'  requires  the  inverse  Fourier  trans- 

2  2       2    2 

forms  of  TSa(coT/2)  and  oj  /([B/2  )   +  w  )  .   The  first  inverse 

transform  has  already  been  determined  and  is 
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t 


/  -  +1,    -T  <  t  <  0 

F-1[TSa^T/2)]   =  (3.59) 

'  -  £+1,  0  <  t  <  T 


2 

Next,  the  inverse  Fourier  transform  of  |H'(u))  |   must  be  deter- 

2 
mined.   Since  |  H  '  (oj  )  |   was  obtained  by  multiplying  H'(u))  by 

-1  i      i  2 
H'(-u)),  F   {|H'(u))|  }  can  be  written  in  the  following  form 


F  1[|H'  (a))  I2]   =   F'1[H'(co)]   *   F"1[H,(-03)]      (3.60) 


Now,  using  the  Fourier  transform  identity  F   (H(-cu)}   <=>  h(-t), 
Equation  3.60  becomes 


F  1[|H'  (a))  |2]   =   [5(t)  -|e"Bt/2y(t)]  *  [6  (-t)  |eBt/2y  (-t)  ] 

(3.61) 


where  the  necessary  inverse  Fourier  transform  has  been  deter- 
mined in  Equation  3.45.   Substituting  this  into  the  convolution 
integral  yields 


F   1[|H'  (go)  |2]       =  /    [fi(t-x)    -|e   B(t   x)y(t-x)] 


[6  (-x)    -|e    B(    x)/2]dx  (3.62) 


and  performing  the  convolution  specified  in  Equation  3.62, 

2 
the  inverse  Fourier  transform  of  |H'  (oi)    is 
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F"1[|H'(aJ)|2]   =   SCt)  -§  e"B|t|/2  (3.63) 


Using  Parseval ' s  theorem,  the  conditional  variance  of 
G'  becomes 

VARfG'ISj^  transmitted]   =   A2TN 

•  [   /  (£+1)  (6(t)  -|eBT/2)dt  +   /  (=£+1)  (6(t)  -|e'Bt/2)dt] 
-T  0 

i  =  0,1  (3.64) 

The  evaluation  of  Equation  3.64  is  fairly  straightforward, 
with  the  exception  that  care  must  be  taken  when  integrating 
the  6 (t)  over  -T  to  zero  and  from  zero  to  T.   There  is  only  one 
singularity  present  at  t  =  0  with  an  area  of  one,  so  that  each 
integral  of  6 (t)  can  be  treated  as  being  equal  to  1/2.   There- 
fore, the  conditional  variance  of  G'  is 


VAR[G'|s.(t)  transmitted] 

=   A2TNQ[^(1  -e"BT/2)]     i  =  0,1  (3.65) 


3 .   Probability  of  Error  Calculation 

With  the  conditional  mean  and  variance  of  G'  specified 
by  Equations  3.53  and  3.64,  it  is  now  possible  to  compute  the 
probability  of  error  performance  of  the  BPSK  coherent  receiver 
with  a  second  order  front-end  filter.   Assuming  the  probability 
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of  the  transmitter  sending  a  "one"  or  a  "zero"  signal  are 
equal  (P(s1(t))  =  P(sQ(t))  =  1/2)),  Equation  2.39  shows  that 
the  average  probability  of  error  P   is 


oo  2 

0'   g 


-mn /v         ~ 

1   9  i      2  .„ 

+    /  -^-  e'X  /Z  dx]  (3.66) 

-oo  /2tt 


where  the  threshold  setting  given  by  Equation  2.33  has  been 
utilized. 

To  compute  P  ,  all  that  is  needed  is  the  specific 
conditional  mean  divided  by  the  square  root  of  the  conditional 
variance  and  using  this  ratio  in  the  limit  of  the  appropriate 
error  function.   First,  for  a  zero  transmitted 

-m       [A2T  -AT/F-  /ITT]  (1  -e"BT/2)2/BT 

— -      =      ^ (3.67) 


g  \/a2TN0(1  -e~BT/2)2/BT 

2 
Recognizing  that  A  T/2  =  E,  the  average  energy  per  bit,  and 

substituting  this  in  Equation  3.67  yields 


—5-   =   /2E/Nn  -  /(2E/N.)  (P  -/E)  (  \/(l  -e  BT/2)2/BT) 
/—  U  U    n: 

g  (3.68) 


Finally 
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-m 


-      =      /2E/Nn(l    -    /P~7e)  (    v/(l  -e    BT/2)J-)  (3.69) 


0N  nj/iJ/  v        v"  'BT 

and   similarly 


-m 


— +1)    v/(l-e-BT/2^ 


=      -    /2E/Nn(/P     ./E  +1)    V(l-e    °i/^)5^  (3.70) 


/v"  °        n^  BT 

g 


The  probability  of  error  performance  can  now  be  determined 
versus  E/Nfl ,  the  signal  to  noise  ratio;  P  . /E ,  the  jammer  to 
signal  ratio;  and  BT,  the  relationship  of  the  filter  bandw.idth 
to  the  inverse  of  the  signal  bandwidth,  by  using  Equations  3.69 
and  3.70  as  the  limits  of  integration  in  Equation  3.66. 

To  evaluate  the  probability  of  error  performance  of 
the  BPSK  coherent  receiver  as  the  filter  bandwidth  approaches 
zero,  L'Hopital's  rule  must  be  utilized.   Applying  this  rule 
to  Equations  3.53  and  3.65  and  then  substituting  the  results 
into  Equation  3.67  shows  that,  as  the  filter  bandwidth  (multi- 
plied by  the  bit  length  T)  goes  to  zero,  the  probability  of 
error  performance  of  the  filtered  coherent  BPSK  receiver 
becomes  that  of  the  coherent  BPSK  receiver  in  the  presence  of 
noise  and  jamming  derived  in  Ref.  2.   Furthermore  when  the 
jammer  to  signal  ratio  goes  to  zero,  the  probability  of  error 
performance  is  the  same  as  that  of  a  coherent  BPSK  receiver  in 
the  presence  of  noise  along  as  derived  in  Ref.  1. 

As  one  more  check,  when  the  filter  bandwidth  approaches 
infinity  both  Equations  3.69  and  3.70  become  zero  and  the 
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probability  of  error  is  one  half.   This  exhibits  similar 
behavior  in  the  limits  of  BT  as  in  the  calculations  performed 
using  the  ideal  filter.   In  a  follow-on  chapter,  the  behavior 
of  the  BPSK  receiver  with  either  type  of  front-end  filter 
will  be  compared  over  the  complete  range  of  variables. 
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IV.       FREQUENCY    SHIFT    KEY    COHERENT    RECEIVER 
WITH    A    FRONT-END    FILTER 

A.   FSK  WITH  A  SINGLE  IDEAL  FILTER 

As  done  with  the  BPSK  receiver,  the  FSK  receiver  will  be 
evaluated  using  an  ideal  front-end  filter,  first  to  gain  in- 
sight into  the  problem  and  second  to  obtain  an  approximation 
to  the  performance  achieved  using  a  higher  order  filter.   In 
the  first  case,  the  front-end  filter  will  be  modelled  by  an 
ideal  bandstop  filter  with  a  center  frequency  ±co   radians/ 
second  and  the  bandstop  region  defined  from  ±co  -B/2  to  ±oj  +  B/2 
where  B  is  the  bandwidth  in  radians/second. 

In  order  to  obtain  a  coherent  FSK  receiver  from  the 
generic  coherent  receiver  of  Figure  A. 2,  s^(t)  is  set  to 


s,(t)       =      -2A  sin  =-((*),  ~wrt)  t  sin  ^(uk  +co~  )  t  (4.1) 

a  2      10  z      l      u 

0    <    t    <    T 


The  input  FSK  signals  are 


s-,  (t)       =      Acosw-.t  0    <_  t    <_  T 


and 


sQ(t)       =      AcostoQt  0    <_  t    <_  T  (4.2) 
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The  optimum  jammer  derived  in  Ref.  2  for  use  against  the 
coherent  FSK  receiver  (without  a  front-end  filter)  is 

2   .1,      4,   .  1 


n  (t)   =   -/P  •  —  sin  T(cu,-u)n)  t  sinT(u).. +W- )  t  (4.3) 

3  n3  /T      2   10       210 

0  <_   t  <  T 

Finally,  the  additive  noise  is  modeled  as  a  white  Gaussian 
noise  process  with  a  power  spectral  density  level  of  Nn/2 
watts/hertz . 

The  first  step  in  determining  the  probability  of  error 
performance  (as  previously  done)  involves  determining  the 
statistics  of  the  random  variable  G'  at  the  output  of  the 
receiver.   Since  it  can  be  shown  that  G1  is  a  conditional 
Gaussian  random  variable,  all  that  is  needed  are  the  conditional 
mean  and  variance  under  the  assumption  that  a  "one"  or  a  "zero" 
are  transmitted.   Finally,  the  probability  of  error  performance 
will  be  calculated  by  applying  these  statistics  to  Equation 
2.39. 

1 .   Calculation  of  the  Conditional  Mean 

In  order  to  determine  the  conditional  mean,  Equation 
2.7  becomes  the  starting  point.   Equation  2.7  is  repeated  here 
for  convenience,  that  is 


E[G|si(t)  transmitted]   =   (s  j,  sd)  +  (n  !  ,  sd)  +^-[  |  |  sQ  |  |  2  -  |  |  s1  |  |  2] 


It  has  already  been  shown  in  the  PSK  system  analysis  that 
for  equal  duration,  equal  amplitude  signals,  the  term 
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2         2 

1/2 ( |  I sQ |    -  |  | s,  |   )  is  zero.   The  frequency  domain  equations 

for  (s?,s-,)  and  (n!,sj  are  given  by  Equations  2.14  and  2.18 
respectively.   Modifying  Equation  2.14  in  order  to  account 
for  the  fact  that  H(co)  is  an  ideal  filter,  (s!,s,)  becomes 


1 
(s!,s,)       =      ~— [      /      S.(-w)S,     (oo)doo    - 

id  2tt         j  i  dp 


dp 

oo    +B/2  -oo    -B/2 

r  o    /      \ ,-.      /    \  j  r  S.(-oj)S,    (oj)dco] 

J  S .  (-oo)  S,     (oo)  doo    -  J  l  dp 

-oo    -B/2  X  Qp  -co    -B/2 

c  c 

i   =    0,1  (4.4) 


In  order  to  calculate  (s!,s,),  Equation  4.4  shows  that  the 

product  of  S.(-co)  and  S,  (oo)  is  needed.   Since  s.(t)  is  defined 
i  dp  i 

over  0  <  t  <  T,  it  can  be  shown  that  S.  (co)  given  by 

*     ,\  ATfc  ,     ,T  -J(w-U)  )T/2 

S  .  (oo)   =   -~-[Sa(co-oo.  )  ~e 
1  Z  1  z. 

-j  (oo+oo  )T/ 2 
+  Sa(oo+aoi)^e       1     ]    i  =  0,1      (4.5) 


S,  (oo)  is  determined  to  be 
dp 

AT  T  -j(a)-oo1)T/2  T  -j(a)+oo1)T/2 

Sd  ^   =   ~2"[Sa  (oo-oj,  )  ^  +Sa  (co+co,)-^ 

T  -j  (oo-ooQ)T/2  T  -j  (oo-oo  )T/2 

-  Sa(co-oo0)^e  -Sa  (go+oOq  )  je  ]       (4.6) 


Before  investigating  the  specific  case  (s*  ,s,)  ,  it 
will  be  convenient  to  look  at  the  relationship  between  the 
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center  frequency  of  the  filter  and  the  signal  and  jammer  fre- 
quencies.  Since  the  filter  contains  a  single  bandstop  region 
and  the  jammer  has  a  spectrum  similar  to  that  of  S,  (co) ,  the 
filter  center  frequency  should  be  located  between  the  jammer 
frequencies  in  order  to  notch  out  the  jammer.   Also,  the 
filter  should  have  a  bandwidth  sufficient  to  cover  both  u)fl 
and  a),  .   Although  the  midpoint  of  coq  and  co,  is  given  by 
(co0+co,  )/2,  it  will  prove  to  be  much  more  convenient  to  define 

co,  and  a),  in  terms  of  co  ,  that  is 
10  c 


CO,     =    CO    +  CO, 

1       c     d 


co„    =   CO    -  CO, 

0       c     d 


(4.7) 


where  co,  =  (co  -con)/2  and  co   =  (co,+co_)/2. 


Substituting  Equation  4.7  for  co,  and  co_  and  multiplying 


S  ,  (co)  by  SA  (-co)  yields 
dp         0 


S   (aOSoC-a)   =  *V 


T  T    jWdT  T        ,  T    "jWdT 

[Sa  (co-co   +co  ,)  7,-Sa  (co-co    +co,)^re  +Sa  (co+co   -co  ,)  Tr-Sa  (co+co    +oo,)^e 

cd2  cd2  ca2  c      a    z 


2  T  2  T 

-Sa     (co-co    +oj,)  rr  -Sa     (co+co    -co,)^-]  (4.8) 

c      d    2  c      d    2 


Before  the  substitution  of  Equation  4.8  into  Equation 

4.4  is  made,  some  simplification  can  be  achieved.   The  contribu- 

2 

tion  of   the   Sa  (x)     functions  with   arguments   of    (co±co    ±oo,)  ,    to 
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the   integrals  with   limits   of    ±u>    ±B/2    is   negligible    for    large 
a)      and   small    B   compared   to   00    .      This    also   holds    true    for   the 
products    Sadu+oj   +(i),)T/2'Sa((j)+o)   -ooJT/2    and   Sa(oo-oo   +00J  T/  2 
•S    (to— oj   -oo,)T/2.      Therefore   substituting  Equation    4.8    into 
Equation   4.4   yields 


,     A2T2      J^dT  °° 

(s0'sd)       =      2~F  ~~ T~ [e  (       f    Sa(a,-wc+wd)^sa(w-ajc-wd)^do 


00    +B/2 

C  rp  rp 

I  Sa  (00-00    +oa  J  ^-Sa  (00-00    -oo,)-~<loo) 

a)   -B/2  c      d    2  c      d    2 

c 


J    d  c  T  T 

+    e  (       /    Sa  (00+00    -03  , )  ^rSa  (00+00    +oo,)7rdoo 

„;  c      d    2  c      d    2 


-00    +B/2 

m  rp 

J  Sa  (00+00    -00 ,)  ~-Sa  (00+00    +ooj^rdoo) 

„  ,J  c      d    2  c      d    2 

-00    -B/2 

c      ' 

00    +B/2 

r         2  T  r  2  T 

-  (      J    Sa  (00-00    +oo,)ydoo    -  J  Sa  (00-00    +oi  J^dw) 

C      a    ^  00    -B/2  C      d 

-00    +B/2 
00  f 

,        2  T  c  2  T 

-  (       /    Sa  (00+00    -go    )ydo3    -  /  Sa  (00+00    -oo,)~doo)]  (4.9) 

c      a    z  -a)  -b/2  C 

c   ' 


Equation  4.9  is  certainly  cumbersome,  but  both  the 
number  of  variables  and  the  complexity  can  be  reduced  by  making 
the  following  change  of  variables  on  the  positive  frequency 
integrals : 
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1)  let   x   =    (o)-o)    )T/2    then    (2/T)dx   =   dw 

2)  when  u)   =    ±°°,    x  =    ±°° 


3)       when  co    =   wc  +  B/2,    x   =    ±BT/4 . 


For    the   negative    frequency    integrals    the   following    substitu- 
tions will   be  made: 

1)  let   x   =    (oo+oo    )T/2,    then    (2/T)dx   =   doo 

2)  when   w   =    ±°° ,    x   =    ±°° 

3)  when  oj  =  -to  ±B/2,  x  =  ±BT/4  . 

After  performing  these  substitutions,  it  is  seen  that  the 
first  two  integrals  are  equal  in  amplitude,  but  have  opposite 
phase  angles.   As  a  result  of  this,  (s',s,)  becomes 


-A2T     r°°   2 
(s0'Sd)   =   ~4ir[(  I    Salx-u>dT/2)dx 


BT/4   2 
/     Sa  (x-a)-,T/2)dx) 
-BT/4 

00      2  BT/4       2 

+    (       /    Sa(x+oo  ,T/2)dx    -  /  Sa    (x-co    T/2)dx) 


—  oo 


-BT/4 


2  cosai.T(       /    Sa  (x+oo  ,T/2)  Sa  (x-oj  ,T/2)  dx 


0         ^  ^         «-X-.wd 


—  OO 


BT/4 
/  Sa(x+u)    T/2)  Sa(x-w   T/2)dx  (4.10 

-BT/4  a  d 


58 


By  similar  analysis 


,2       °°  BT/4 

(s'  s  )   =  -j£-[  (   /  Sa(x-oj-,T/2)dx  -      /     Sa2(x-oo  T/2)  dx! 

-~         a         -BT/4  a 


00   2  BT/4   2 

+  (   /  Sa(x-co  T/2)dx  -      /    Sa  (x+oj  ,T/2)  dx) 
a  -BT/4  a 


-2cosa)0T(   /  Sa(x+wdT/2)Sa(x-GodT/2)dx 

—00 

BT/4 
/      Sa  (x+oo,  T/2  )Sa(x-oo,  T/2 )  dx)  ]  (4.11) 

-BT/4  a  a 


The  next  term  in  the  conditional  mean  is  (n!,s,) 

D   a 

which  is  given  by  Equation  2.18.   As  previously  done,  this 
equation  is  modified  to  become 


1 
(n!,  s,)   =  ~[   /  S  ,  (oj)N.  (-oj)doj  • 
:   d      2tt  _J       dp     : 

ooc+B/2 

/       S,  (oj)N.  (-co)daj 

co  -B/2         dP     ^ 
c 

-co  +B/2 
c   / 

/        S,  (co)N.  (-co)doo]  (4.12) 

-co  -B/2'  dP     ^ 

c 


Since  H(co)  is  an  ideal  filter,  in  order  to  form  the  product 

of  N.  (-co)  and  S,  (oo)  ,  the  spectrum  of  n.  (t)  must  be  found. 
Z)  dp  j 

The  equation  for  n.(t)  (see  Equation  4.3)  shows  a  truncated 
process  similar  to  that  found  for  the  BPSK  case.   Thus  N.  (co) 
is  determined  to  be 
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m  m  "j  (W-O),  )T/2 

N.  (u>)   =   /Pnj/T  j[Sa(oo-a)1)|e 

T  -j  (w+(Jo1)T/2 
+  Satw+co-^^e 

-j(o)-oj  )T/2  -j(oo-co  )T/2 

-  Sa(oo-oo0)^e       u      -  Sa  (co+u)  )  |e      u     ]  (4.13) 


Now,  multiplying  N.  (-w)  by  S,  (ou)  ,  where  S,  (00)  is  given  by 
Equation  4.9  yields  eight  terms.   Four  are  on  the  positive 
side  of  the  co  =  0  axis  and  four  are  on  the  negative  side  of 
the  w  =  0  axis.   For  convenience  just  the  positive  frequency 
terms  are  shown  here 


2 
N.(-a>)S d  (u>)   =   ^h  /Pnj/T[Sa2(o3-o)0)| 


T         T     2       T 
-  2  cos  todTSa  (oo-co_)  jSa  (cu+oo,  )  y  +  Sa "w-io,  )  y]  (4.14) 


Expressing  con  and  u),  as  co  ±00-,  as  before,  and  substituting 
Equation  4.14  into  Equation  4.12  under  the  assumptions  stated 
in  the  development  of  (s!,s,),  yields  an  expression  similar 
to  Equation  4.9.   The  difference  here  is  that  the  Sa(co-w  +co,)T/ 

2Sa(a3-oo  -io,)T/2  terms  are  multiplied  by  2  cos  u),  T  rather 

j5dT 
than  e     .   By  making  the  same  change  of  variables  as  in 

Equation  4.10  and  recognizing  that  Equation  4.14  covers  only 

half  of  the  frequency  spectrum  (which  is  symmetrical) , 

(n ! ,s  ,)  becomes 
3      a 


60 


GO 

ni'Sd}       =      27AT/Pn i/T[(      /    Sa2(x-wdT/2)dx 

■J  -J  _0O 


BT/4 

/  Sa^(x-(jo,T/2)dx)  (      /    Sa^  (x+oo  ,T/2) 

-BT/4  a  -oo  a 

BT/4       2 
/  Sa    (x+w    T/2)dx) 

-BT/4  a 


-2coscodT(      /      Sa(x+oodT/2)  Sa(x-wdT/2)dx 


— oo 


BT/4 
/  Sa(x+u),T/2)Sa(x-a),T/2)dx]  (4.15) 

-BT/4  a  a 


12  2 

Finally,  with  (s|,sd),  (n£,sd)  and  ^" (  j  I  sQ  j  |   -  |  |  s1  |  |  \ 

determined,  the  conditional  mean  of  G'  can  be  expressed  as 


2 

E[G'|s.    transmitted]       =       [AT/P    ./T   +-l"1+1   K^]  -^- 
1     l  nj  2      2tt 

oo  BT/4 

•[(      /    Sa2  (x+o)    T/2)dx   -  /  Sa2  (x+ca    T/2 )  dx) 

a  -BT/4 


°°      2  BT/4       „ 

+    (       /   Sa    (x-oo.T/2)dx   -  /  Sa     (x-to,T/2)  dx) 

Q  -BT/4  a 


-2cosoo,T(       /      Sa  (x-oo  ,T/2)  Sa  (x+oo  ,T/2)  dx 


BT/4 
/  Sa  (x-o)   T/2)Sa(x+co   T/2)dx)  ]  ,       i    =    0,1  (4.16) 

-BT/4  d  d 
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2.   Calculation  of  the  Conditional  Variance 

The  conditional  variance  of  G'  in  the  frequency  domain 
is  given  by  Equation  2.27,  which  can  be  modified  to  take  on 
the  following  form 


N0 
VAR[G'  s.(t)  transmitted]   =  ~^- 

1  <£TT 


to    +B/2 
c 

[    /    |S,    (to)|dco    -  /  |s,    (a))|doo],    i   =    0,1         (4.17) 

0  dP  a)    -B/2  dp 


c 


due  to  the  fact  that  H  (to)  is  assumed  to  be  an  ideal  filter. 

Therefore  in  order  to  calculate  the  conditional  variance,  all 

2 
that  is  needed  is  the  determination  of  |  S  ,  (co)  |  .   Since 

2 

|  S,     (to)  I       =    S,     (co)S,     (-co)     and    S.     (oo)     is    given   by   Equation    4.8, 

the   terms   on   the   right   hand   side   of    the  to    =    0    axis    are    given  by 

2    2 

|  S  .,    (oj)  J         =      — -A — [Sa  (co-co    +oj  J  ~n-Sa  (co-co   -co,)-=- 
1    dp  '  4  c      d    z  c      d    2 

positive 

frequency  -2cosco  ,TSa  (co-to    +oo,)  ^Sa  (co-co    -co,)^]  (4.18) 

terms  d  c      d    2  c      d    2 


where   co.    =   co    -co,    and   co,     =   co    +co,. 
0  c      d  1  c      d 

Since   the   limits    of   the    integrals    in   Equation  4.17    are 

on   the   right  hand   side   of   the    co   =    0    axis,    for   large   co      and 

2 

small  B  compared  to  co  ,  the  terms  of  |  S ,  (co)  |   on  the  left  hand 

side  of  the  co  =  0  axis  can  be  neglected.   By  substituting 
Equation  4.18  into  4.17  it  is  quickly  recognized  that  the  inte- 
grals are  of  the  same  form  as  Equations  4.9,  4.10  and  4.14. 
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So,  instead  of  repeating  the  steps  which  have  been  demonstrated 
before,  the  conditional  variance  is  shown  here  without  further 
development 


2  °° 

VAR[G'|s.(t)    transmitted]       =      —-?-[  (      /    Sa^x-oo  ,T/2)  dx 

1  2tt  j  d 


BT/4  oo 

/  Salx-oo,T/2)dx)    +    (      /    Salx+a)  -,T/2)  dx 

-BT/4  a  a 


BT/4 
/  Salx+oo    T/2)dx) 

-BT/4  a 


2coscodT(      /    Sa(x-oodT/2)Sa  (x+a)dT/2)  dx 


—  oo 


BT/4 
/  Sa(x-oo   T/2)  Sa(x+w,T/2)dx)  ]  ,       i   =    0,1  (4.19) 

-BT/4  a  a 


3.   Probability  of  Error  Calculation 

With  the  conditional  mean  and  variance  of  G'  given  by 
Equations  4.16  and  4.19,  it  is  now  possible  to  compute  the 
probability  of  error  performance  of  the  BFSK  coherent  receiver 
with  an  ideal  front-end  filter.   Assuming  that  the  probability 
of  the  transmitter  sending  a  "one"  or  a  "zero"  signal  are 
equal  (P(s  (t) )  =  P(s  (t))  =  1/2)),  Equation  2.39  shows  that 
the  average  probability  of  error  P   is 


0  -ra,//v~  „ 

oo  2  1         cr  2 

Pe      =      \[  _/    J-  e"X   /2dx   +         /  J-  e"X   /2dx]        (4.20 

-m_//v         /2t\  -°°  /2tt 

0         g 
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where  the  threshold  setting  given  by  Equation  2.3  3  has  been 
utilized. 

In  order  to  compute  P  ,  all  that  is  needed  is  the 
specific  conditional  mean  divided  by  the  square  root  of  the 
conditional  variance  and  using  this  ratio  in  the  limit  of  the 

appropriate  error  function  integral.   Recognizing  that 

2 
A  T/2  =  E,  the  average  energy  per  bit,  and  making  the  appro- 
priate cancellations/  mn//v~~  is  determined  to  be 

-mQ 


=   /E/Nn(l  -  /2Pn./E)  /Filter  Factor         (4.21) 
/v~         U         : 

g 


and  m,  //v-  becomes 
1   g 


"ml 


=   -/E/N  (1 ■+  /2P  ./E)  /Filter  Factor        (4.22) 
/v~  U        nD 

g 

where  the  Filter  Factor  is  defined  as 

Filter  Factor   =  -^~ 

2tt 

°°  9  BT/4 

•  [(   /  Sa  (x-co -T/2)dx  -       /     Sa^x-w  .T/2 )  dx) 
a  -BT/4  a 

°°  9  BT/4 

+  (   /  Sa^(x-w  T/2)dx  -       /     Sa^(x-oo.T/2 )  dx) 
a  -BT/4 

oo 

-    2  cosw    T(       /    Sa  (x-codT/2)  Sa  (x+oodT/2)dx 

—  oo 

BT/4 
/  Sa  (x-co -.T/2)  Sa  (x+w    T/2)  dx)  ]  (4.23) 

-BT/4  a  d 
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The  probability  of  error  performance  can  now  be  determined 
versus  E/NQ ,  the  signal  to  noise  ratio;  P  ,/E,    the  jammer  to 
signal  ratio;  and  BT,  the  relationship  of  the  filter  bandwidth 
to  the  inverse  of  the  signal  bandwidth,  by  using  Equations 
4.21  and  4.22  as  the  limits  of  integration  in  Equation  4.20. 
It  is  desirable  to  compare  these  results  with  those 
previously  obtained  for  both  the  jamming  problem  and  for 
the  general  FSK  coherent  receiver.   In  order  to  compare  the 
results  for  the  probability  of  error  performance  of  the  coher- 
ent FSK  receiver  with  no  filtering,  it  is  necessary  to  analyze 
how  those  results  were  obtained.   For  the  coherent  receiver 
operating  in  just  additive  white  Gaussian  noise  alone,  the 
probability  of  error  from  F.ef.  1  is 


00   1      2/9 
P    =        f   —  e~U  7   du  (4.24 

/2tt 


/E/Nq 


2 
This  result  is  obtained  when  the  energy  per  bit  is  A  T/2 

and  the  normalized  cross  correlation  p"  is  zero,  where  (see 

Ref.  1) 


=      \        J         [yQ(t)    +   yjttjjdt  (4.25) 

fc0 

=      |        /         y0(t)Yl(t)dt  (4.26) 

t0 
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and  tQ ,  t_  are  replaced  by  0,  T  respectively.   From  these 

expressions  it  is  easy  to  show  that  (jo,T  must  equal  m\   and 

2 
00-.T  =  kTr  (where  n  and  k  are  integers)  for  E  =  A  T/2.   For  the 

FSK  modulation  p  is 


p   =   Sa(io  +wQ)T  +  Sa(u),  -coQ)T  (4.27) 


For  large  (oa,+con)T  the  Sa(oo,+oo  )T  is  negligible  and  p  is 
zero  when  (oo,-u)0)T  is  equal  to  miT  (where  m  is  an  integer)  . 
When  this  is  substituted  into  Equation  4.23,  recalling  that 
co-jT  =  (oo, -u)n)  T/2 ,  it  can  be  seen  that  when  BT  is  zero  the 
Filter  Factor  of  Equation  4.23  is  1.0.   For  values  of  oo,T 
that  are  odd  multiplies  of  tt/2  the  Filter  Factor  is  exactly 
equal  to  1.0,  but  for  even  multiples  it  becomes  1  plus  a  small 
quantity  that  is  negligible.   So  again,  as  in  the  PSK  ideal 
front-end  filter  case,  for  BT  equal  to  zero,  Equations  4.53 
and  4.54  are  in  agreement  with  the  results  of  Ref .  2  for  the 
coherent  FSK  receiver  operating  against  the  optimum  jammer. 
For  P  ./E  equal  to  zero  these  results  correspond  to  the  results 
of  Ref.  1  for  the  coherent  FSK  receiver  operating  in  additive 
white  Gaussian  noise  alone. 

B.   FSK  WITH  A  SECOND  ORDER  FRONT-END  FILTER 

In  order  to  evaluate  the  effect  of  a  second  order  band- 
stop  filter  on  the  probability  of  error  performance  of  a 
coherent  FSK  receiver,  the  filter  specified  by  Equation  3.37 
will  be  utilized.   The  use  of  the  equivalent  baseband  filter, 
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along  with  Parseval's  theorem  will  not  greatly  aid  in  the 
determination  of  the  conditional  mean  and  variance  of  G* 
(the  output  of  the  receiver  of  Figure  A. 2).   For  convenience, 
the  parameters  of  the  coherent  FSK  receiver  model  are  repeated 
here  (see  Equations  4.1  and  4.2) 


s2  (t)   =   -2A  sin  -^(co, -oj_)  t  sin  2"(oo,+ooQ)t 


sQ  (t)   =   A  cos  coQt 


s,  (t)   =  A  cos  out      0  <_   t  <_   T 


The  optimum  jammer  derived  in  Ref.  2  is 


n.(t)   =   -/P  .  —sin  ^■(a)1-ujA)  t  sin  ~-(go,  +a)n)  t    0  <  t  <  T 
3  nj7T210        211)        —    — 


and   the   filter   of   Equation    3.37    is 


2  2 

0)        -     0) 
H   0)         = 


2         2 

oo  n  -  GO      +  ]ooB 


where  go   will  be  replaced  by  go   as  in  the  FSK  calculations 
0  c 

involving  the  ideal  filter  case. 

1 .   Calculation  of  the  Conditional  Mean 

As  before,  Equation  2.7  will  be  used  to  determine  the 
conditional  mean  of  G'.   Since  the  FSK  signals  are  equal 
energy  signals,  all  that  is  needed  from  Equation  2.7  are  the 
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terms  (n!,sd)  and  (s!,s,)  with  i  =  0,1.   These  terms  are  given 
by  Equations  2.14  and  2.18.   First  (s!,sd)  will  be  calculated 
by  evaluating  the  product  of  S.  (-co)  with  S,  (oo)  along  with 
H(-co)  and  then  substituting  into  Equation  2.18  which  is 
repeated  here  for  convenience 


oo 

(si'sd}   =   27[  I   S.(-u)Sd  (u)H(-a»)du    i  =  0,1 

-iYt  * 


Investigating  the  particular  case  of  (s',s,),  the  product 
of  S-(-co)  and  S,  (co)  is  given  by  Equation  4.8  where  oo  +co,  and 
oo  -co,  were  substituted  for  co,  and  oo   respectively.   Thus, 
substituting  Equation  4.8  into  Equation  2.18  yields 


2    2 
,     ,  x  1    AT 

lS0,Sd;  2tt        4 


jco,T  °° 

[e  (      /    sa(co-oo   +oo    )TSa(co-co    -co ,  )&  (-co)dco) 

c      d    /  c      a    l 


_ja3dT    r  t  t 

e  (       J    Sa(w+a)    -to    )-=-Sa(co+co    +co    )— H(-oo)doo) 

c      a.    z  c      q   z 


—  oo 


,2  T  2  T 

J    Sa(co-co    +oo -.)^-H(-co)dco    -         J    Sa  (to+co    -to  ,)-=-H  (-co)  dec]  (4.28) 

'  c      d    2  '  c      d    2 


to  further  determine  (s',s, )  it  is  necessary  to  investigate 
both  H(-co  -2x/T)  and  H(oo  -2x/T)  .   By  substituting  co  -2X/T  for 
oo  in  Equation  3.37,  expanding,  and  multiplying  numerator  and 
denominator  by  T/4co   results  in 
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H(o)    -2x/T)       = 
c 

2 

X 

X      a)    T 
c 

(x 

x2        jBT 

jBx 

2co 
c 

(4.30) 


The  filter  transfer  function  was  derived  under  the  assumption 
that  oa   >>  B  so  H  (go  -2x/T)  becomes 


2 
x 
X  - 


GO  T 

H(oo  -2x/T)   =   =-°- (4.31) 

^X  +oo  T;+  4 


and  similarly 


2 
x  + 


0)  T 
H(-u)  -2x/T)   =   — 5-£ (4.32) 

lX   co  T;+  4 
c 


For  ease  of  notation  let  H,  (x)  =  H(co  -2x/T)  and  H„  (x)  =  H(-co  -2x/T! 
Then  separating  H. (x)  and  H~ (x)  into  real  and  imaginary  parts 
yields 


2  0 

i  x  ,  2 

Re{H,(x)}   =   5 - (4.33) 

x   2    BT  2 


"BT-     x  . 

T-(X-STT) 

Im{Hn(x)}   =   =- - (4.34) 

1  ,     x2  2    BT.2 

(X-^T}   +(^ 
c 
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2 

Re{H  (x) }   =   = (4.35) 

c 


-BT  .       ,  x2 , 

Im{H  (x)  }   =   . - (4.36) 

,   J  x  ,  2   ,BT.2 

c 


With  Equations  4.33  through  4.36  it  is  now  possible  to  expand 
Equation  4.29  into  real  and  imaginary  parts  which  yields 


2  oo 

(s£,sd)   =   ^-[cosu)dT[   /  Sa(x-0JdT/2)Sa(x+o)dT/2) 


[Red^  (x) }+Re{H2 (x) }dx 


+    j         /    Sa(x-todT/2)Sa(x+ajdT/2)  (Im{H,  (x)}fIm{H2  (x)  })dx] 


+    sino)dT[      /    Sa(x-o)dT/2)Sa(x+wdT/2)  (Im{H1(x)  }-Im{H2(x)  })dx 


+    j         /    Sa(x-codT/2)Sa(x+codT/2)  (RelHj^  (x)  }-Re{H2(x)  })dx] 


-     [      /    Sa(x+u>dT/2)Re{H    (x)  }dx  +       /    Sa2(x-oodT/2  )  Re{H2  (x)  }dx 


+    j         /    Sa  (x+codT/2)  ImlHj^  (x)  }+Sa2(x-a3dT/2)  Im{H2  (x)  }dx]  ]  (4.37 
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With  the  real  and  imaginary  parts  of  H, (x)  and  H~  (x)  as  given 
by  Equations  4.33  through  4.36,  it  is  easy  to  see  that  the 
imaginary  terms  integrate  to  zero  because  of  the  odd  symmetry 
about  x  =  0  that  the  functions  of  Equation  4.37  exhibit.   Thus 
(Sq,s,)  becomes 


(s0'Sd}   =   "TfT"[   /  Sa2(x+codT/2)Re{H1(x)  }dx 


+    /  Sa2(x+co,T/2)Re{H  (x)}dx 


-  coso)dT    /  Sa  (x-oodT/2)Sa  (x+o)dT/2)  (Re{H1  (x)  }+Re{H2(x)  })dx 

— oo 

oo 

-  sinco    T        /    Sa  (x-o)dT/2)Sa(x+o)dT/2)  (Im{H1(x)  }-Im{H2(x)  })dx] 

—  00 

(4.38) 
and  by  similar  analysis 


(si'sd}       =      TF[      /    Sa2(x+wdT/2)Re{H1(x)  }dx 


°°      2 
+         /    Salx-ojdT/2)  Re{H2  (x)  }dx 

— oo 

oo 

-  cosojdT         /    Sa  (x-ojdT/2)  Sa  (x+oadT/2)  (Re{Hx  (x)  }+    Re{H2(x)})dx 

—  oo 

oo 

-  sincodT         /    Sa  (x-oodT/2)Sa  (x+oodT/2)  (ImlHj^Cx)  }-Im{H2  (x)  }dx] 

—  oo 

(4.39) 
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The  next  term  that  must  be  computed  is  (n!  ,s,),  given 
by  Equation  2.18  and  repeated  here  for  convenience 


(n!,s,)   =  i-   /  S ,  (u>)N.  (-u))H(-u>)da> 
j   d      2tt    ;   dp    .  j 


—  00 


As  before,  the  calculation  of  this  term  starts  with  the  pro- 
duct of  N.(-co)  and  S,  (co)  .   This  product  when  fully  expanded 
contains  sixteen  terms.   Eight  are  products  of  Sa (x)  functions 
centered  on  opposite  sides  of  the  co  =  0  axis  and,  for  large  u), 
oon ,  are  essentially  equal  to  zero.   By  expressing  con  as  to  -co, 
and  co-  as  co  +co , ,  substituting  into  Equation  2.18,  and  performing 
the  same  change  of  variables  as  in  the  expansion  of  (s!,s,), 

it  can  be  shown  that  (n!,s,)  becomes 

j      d 


CO 

(n!,s,)       =      ^/P     ./T[       /    Sa2(x+co,T/2)  (H(co    -2x/T)+H(-co    -2x/T))dx 
id  4tt      nn  J  d  c  c 

->  •>  — OO 


°°      2 
+         /    Salx-co    T/2)  [H(co    -2x/T)+H(-co    -2x/T))dx 


00  "J^^T 

-    2  f    Sa  (x-coJT/2)Sa  (x+co -.T/2)  [e  H  (co    -2x/T) 


jajdT 

+    e  H(-coc-2x/T)  )dx]  (4.40) 


The  integrals  in  Equation  4.40  are  now  of  the  same 
form  as  those  in  Equation  4.29.   Therefore,  the  same  change  of 
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variables  operation  can  be  performed  as  done  in  the  derivation 
of  Equation  4.40.   This  results  in 


(n!,sd)   =   ||/P./T[   /  Sa2(x-wdT/2)  (Re{H1(x)  }+Re{H2(x)  })dx 

■J  J      _oo 


+    /  Salx+o)dT/2)  (Re{H,  (x)  }+Re{H  (x)  } )  dx 

—  CD 

oo 

-  2coso)dT    /  Sa(x-oa  T/2)Sa(x+a)  T/2)  (Re{H1  (x)  }+Re{H2  (x)  })dx 

— oo 

oo 

-  2sino)dT   /  Sa(x-a)  T/2)Sa(x+o)dT/2)  (Im{H1(x)  }-Im{H2(x)  }dx] 

-co 

(4.41) 

where  Re{H, (x)},  Re{H2(x)},  Im{H,  (x)  }  and  Im{H2(x)}  are  defined 
by  Equations  4.3  3  to  4.36.   The  conditional  mean  of  G*  under 
the  assumption  that  sn(t)  or  s->  (t)  was  transmitted  can  be 
expressed  as  the  sum  of  Equation  4.38  or  4.39  with  Equation 
4.41. 

2 .   Calculation  of  the  Conditional  Variance 

The  conditional  variance  of  G*  is  given  by  Equation 
2.27  and  repeated  here  for  convenience 


N      °° 
VAR[G'|si(t)  transmitted]   =   ~    /  |  Sd  (oo)  |  2  |  H  (co)  |  dou 


2 

by  recognizing  that  j  H  ( oo )  |   is  H  (co)  H  (-co)  ,  for  real  h(t), 

Equation  2.27  can  be  expressed  as 
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N      °° 
VAR[G,|si(t)  transmitted]   =   47   /  I  Sd  (w)  |  2H  (00)  H  (-00)  dw  (4.42) 


—00 


In  order  to  begin  the  calculation  of  the  conditional  variance, 

1       1  2 
as  for  the  ideal  filter  case,   S,  (co)    must  be  determined  by 

1  dp    '  2 

multiplying  S,  (to)  by  S,  (-00)  .   This  product  contains  sixteen 

terms  as  encountered  in  the  product  of  N.(-co)  and  S  ,  (co)  . 

j  dp 

2 

Fortunately,  |  S ,  (go)  |   only  differs  from  N.(-co)S,  (co)  by  a 

constant.   Therefore  it  is  possible  to  perform  a  similar  set 
of  operations  as  done  in  the  derivation  of  Equation  4.29  so 
as  to  arrive  at 

NQA2T 

VAR[G'|s.  transmitted]   =   — -. 

1     1  4tt 


CO 

[       /    Sa^x-co.T/2)  [H  (2x/T-co    )H(co    -2x/T)+H(-2x/T-co    )H(uj    +2x/T)  )dx 


00 

+    [       /    Sa^x+co    T/2)  (H(2x/T-co    )H(oo    -2x/T) +H  (00    +  2x/T)  H  (-00    -2x/T)  )  dx 


00  j^n"1" 

-2         /    Sa(x+todT/2)  Sa(x-codT/2)  (e  H  (  2x/T-a>c)  H  (ooc"2x/T) 


+    e         a   H(2x/T+oj    )H(-2x/T-co    ))dx]  i    =    0,1  (4.43) 


where  H(-co  -2x/T)  and  H  (co  -2x/T)  are  given  by  Equations  4.31 

and  4.32.   It  can  be  demonstrated  that  H(2x/T-co  )  and 

c 

H(2x/T+co    )    are   given  by 
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(x  +■ 

2 
X    ) 
c 

(x 

2 

+    *«, 

to    T 
c 

)_  JBT 
;       4 

H(2x/T+ooc)   =         q  (4.44) 


and 


2 

/    x  . 
(x ^) 

H(2x/T-co  )   =   » (4.45) 

(x-  *)_JBT 

c 


By  substituting  Equations  4.30,  4.31,  4.44  and  4.45  into 
Equation  4.42  and  expanding  it  into  real  and  imaginary  parts 
the  conditional  variance  of  G'  becomes 

N  A2T 

VAR[G'|s.    transmitted]       =      — —. 

'1  4tt 


2  x2     2 

(x-JTt)  (x+STt) 


*-\  UU  J- 

/    Salx-co    T/2)  [ = + 


c 


x2     2       ,BT\2  ,  x2 . 2       ,BT.2 

(x-^Tt)     +(T)  (x+^Pr)     +(T} 

c  c 


x2     2  x2     2 

w,t         (x  "[^T)  (x  +^rT} 

+      /    Sa(x+-^-)  [ ^— 2 +    ^ ]dx 

,  x    ,2       ,BT\2  ,  x     >2       ,BT.2 

(X"^       +  {~T]  {X+^]      +  (TT) 

c  c 

x2     2  x2     2 

0JdT  oo    T  (X"VT)  (X  +  ^TT} 

-2cosoj    T         /    Sa(x — ^— )  S     (x+-=— ) [ =— + 5 ]  dx 

a      -oo  Z  .  x     ,  2        ,BT\2       ,  x     ,  2        ,BT>2 

^"jtt5    +  (-t}      (x+^rT}    +(^ 

c  c 
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x2     2  2     2 

o)dT         oo.t       (x_^nr)  {x+zr~r) 

-2cosco,T      /    Sa(x — 2_)  Sa  (x+-S-)  [ 5— + j— ]dx} 

-«  ,        x     ,  2     .BT.2       ,        x     , 2^,BT,2 

c  c 

i    =    0,1  (4.47) 

3 .   Probability  of  Error  Calculation 

The  conditional  mean  of  G1  is  given  by  the  sum  of 
Equation  4.41  with  either  Equation  4.38  or  4.39  depending 
whether  sfi(t)  or  s, (t)  were  transmitted,  respectively.   The 
conditional  variance  is  Equation  4.47.   With  the  conditional 
mean  and  variance  determined  it  is  now  possible  to  compute  the 
probability  of  error  performance  of  the  FSK  coherent  receiver 
with  an  ideal  front-end  filter.   Assuming  the  probability  of 
the  transmitter  sending  a  "one"  or  a  "zero"  signal  are  equal 
(P(s1(t))  =  P(sQ(t))  =  1/2)),  Equation  2.39  shows  that  the 
average  probability  of  error  P   is 


2„        "V1^   .      2/0 

Pe  -  il     ^J     -L-   e"x  /2dx  +    /        A-e^  /2dx] 
-mn//v~        /2~7T  -°°  /2~¥ 

0    g  (4.48) 


where  the  threshold  setting  given  by  Equation  2.33  has  been 
utilized. 

In  order  to  compute  P  ,  all  that  is  needed  is  the 
specific  conditional  mean  divided  by  the  square  root  of  the 
conditional  variance  since  this  ratio  must  be  used  in  the  limit 
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of  the  appropriate  error  function.   However,  the  integrals  of 
Equation  4.4  7  do  not  match  those  of  the  conditional  mean,  as 
found  in  every  case  investigated  thus  far.   So  before  going 
on  to  compute  the  limits  of  the  error  functions  of  Equation 
4.48,  the  behavior  of  the  conditional  means  and  the  variances 
will  be  investigated  with  respect  to  id  T.   Referring  to  Equa- 
tions 4.33  through  4.36  and  4.44  to  4.45,  each  has  a  term  of 

2 

the  form  x  /to  T  in  both  the  numerator  and  denominator. 

Remembering  that  the  filter  was  derived  using  an  assumption 

2 

of  large  go  ,  the  effect  of  the  x  /oj  T  terms  can  become  negligi- 
ble when  used  in  the  integrals  involving  products  of  two  Sa (x) 
functions,  since  these  functions  decrease  rapidly  with  increas- 
ing x.   Therefore  for  large  oo  T,  a  Filter  Factor  can  be  defined 
as 


°°      oj,T        0JdT        2 
Filter  Factor   =     /  Sa  (x+-^-)  +Sa  (x — ^-)  ]  (  2  XBT  2)dx 

-«  x    +(-2") 

gj,T  oo,T  2 

-    2cosoadT         /    Sa(x — 2_)Sa(x+-2-)  (^ X£T    2)dx  (4.49) 

x  +  C-^-) 


Thus  the  conditional  variance  becomes 

A2TNQ 

Var[G'|s.  transmitted]   =   — ^ —  x  Filter  Factor  (4.50) 

1  l  2tt 

i  =  0,1 


and  the  conditional  mean  becomes 
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E[G'|s.  transmitted]   =   [AT/P  ./T  -11+1  ^-] 
1  1  ny  2 

xJLxFilter  Factor    i  =  0,1      (4.51) 

The  Filter  Factor  was  computed  via  numerical  integration 
methods  and  compared  with  the  complete  expressions  of  Equation 
4.47  for  u)  T  as  small  as  50  and  as  large  as  1000  with  negligi- 
ble errors.   An  example  of  a  large  u  T  would  be  a  communications 
system  operating  at  1.0  MHz  with  a  bit  length  of  13.3  milli- 
seconds (75  baud)  .   The  oo  T  term  would  have  a  value  of  13,300 

2 
which  further  supports  dropping  the  x  /w  T  term  from  the  equations. 

Evaluating  now  mn//v  by  dividing  Equation  4.51  by 

the  square  root  of  Equation  4.50  where  the  Filter  Factor  is 

defined  by  Equation  4.49,  iru/Zv-  becomes 


mn       [AT/P  ./T  -  A  T/2]  x^-  x Filter  Factor 

-°-      =      n=>  2tt (4.52) 

vg  va  TN„x^-  xFilter    Factor 

0  2tt 

2 
Recognizing  that  A  T/2  =  E,  the  average  energy  per  bit,  and 

substituting  this  in  Equation  4.52  yields 


-mQ 


=   /E/Nn(l  -  /2P  ./E  /1/2-tt  x  Filter  Factor         (4.53) 


0V        nj 

/v  J 
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and  similarly 


"m!        

— -  =      -    /E/Nn(l  +  /2P  ./E)  /1/2tt  x  Filter  Factor      (4.54) 
/ —  0         ni 
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The  probability  of  error  performance  can  now  be  determined 
versus  E/NQ,  the  signal  to  noise  ratio;  P  ■ /E ,    the  jammer 
to  signal  ratio;  and  BT  (the  relationship  of  the  filter  band- 
width to  the  inverse  of  the  signal  bandwidth) ,  by  using 
Equations  4.53  and  4.54  as  the  limits  of  integration  in 
Equation  4.48. 

As  with  the  ideal  filter  case,  in  order  to  compare  the 
results  for  the  probability  of  error  performance  of  the 
coherent  FSK  receiver  with  no  filtering,  it  is  necessary  to 

remember  that  those  results  were  obtained  using  p  =  0  and 

2  — 

E  =  A  T/2,  where  p  and  E  are  defined  by  Equations  4.2  6  and 

4.25  respectively.   Therefore,  for  large  (oo,+oon)T,  the 

Sa(oo, +wn)  T  term  is  negligible  and  p  is  zero  when  (co.-io-JT  is 

equal  to  m  .   When  this  is  substituted  in  Equation  4.23, 

recalling  that  co,T  =  (w,-to_)T/2f  it  can  be  seen  that  when  BT 

is  zero  the  Filter  Factor  is  2tt.   For  values  of  (jo-,T  that  are 

d 

odd  multiples  of  tt/2  the  Filter  Factor  is  exactly  equal  to 
2tt,  but  for  even  multiples  it  is  equal  to  2  plus  a  quantity 
which  is  negligible.   So  again,  as  in  the  FSK  ideal  front-end 
filter  case,  for  BT  equal  to  zero,  Equations  4.53  and  4.54  are 
in  agreement  with  the  results  of  Ref .  2  for  the  coherent  FSK 
receiver  operating  against  the  optimum  jammer.   For  P  ./E 
equal  to  zero,  these  results  correspond  to  the  results  of  Ref. 
1  for  the  coherent  FSK  receiver  operating  in  additive  white 
Gaussian  noise  alone.   Again,  as  with  the  coherent  BPSK  receiver, 
as  the  BT  approaches  zero,  the  behavior  of  the  probability  of 
error  performance  of  the  coherent  FSK  receiver  calculated  using 
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a  second  order  front-end  filter  is  similar  to  the  performance 
calculated  using  the  ideal  front-end  filter.   Later  the  per- 
formance comparison  will  be  performed  throughout  the  range  of 
pertinent  values  for  the  variable  parameters. 
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V.   RESULTS  AND  DISCUSSION 

A.   GENERAL 

For  each  of  the  four  receiver  configurations  (PSK  and  FSK 
receivers  with  both  ideal  and  second  order  front-end  filters) , 
the  results  for  the  probability  of  error  performance  have  been 
calculated  and  plotted  as  functions  of  P  ./E    (jammer  to  signal 
ratio) ,  E/NQ  (signal  to  noise  ratio)  and  BT  (the  filter  band- 
width, bit  duration  product).   The  meaning  of  the  jammer  to 
signal  and  signal  to  noise  ratios  are  self  explanatory. 
However,  the  term  BT  deserves  some  attention.   Since  T  is  the 
time  duration  of  s,  (t)  and  sn(t)  ,  the  spectra  of  each  have 
been  shown  to  be  functions  of  the  form  Sa(coT/2)  where  Sa(x) 
has  been  defined  as  sin(x)/x.   The  Sa(ooT/2)  function  has  nulls 
when  the  value  of  ooT/2  is  an  integer  multiple  of  n   or  when 
oo  =  N(2tt/T)  . 

From  the  Sa(ooT/2)  spectrum  it  can  be  shown  that  the  value 
of  B  =  2tt/T  represents  approximately  the  3  dB  (or  half  power) 
bandwidth  of  the  signal  (actually  2tt/T  is  the  4  dB  bandwidth)  . 
The  term  BT  can  thus  be  thought  of  as  the  ratio  of  the  filter 
bandwidth  to  the  signal  bandwidth  for  intuitive  purposes.  It 
is  still  convenient  to  refer  in  further  discussions  to  the 
location  of  the  nulls  of  the  Sa(ojT/2)  function. 

The  analysis  of  the  calculated  probability  of  error  per- 
formance for  each  receiver  configuration  uses  three  types  of 
plots.   The  first  is  the  familiar  plot  of  probability  of 
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error  versus  signal  to  noise  ratio,  with  constant  values  of 

BT  and  P  ./E.   The  values  of  P  . /E  used  are  0,  .1,  .5,  1,  and 

nj  r\y  '         '         '       ' 

10.   These  values  correspond  to  no  jammer  present,  and  jammer 
to  noise  ratios  of  -10  dB,  -3  dB,  0  dB,  and  10  dB .   The  par- 
ticular values  of  BT  depend  on  the  modulation  structure 
analyzed,  and  will  be  discussed  later.   The  second  type  of 
plot  will  allow  the  analysis  of  the  probability  of  error 
performance  to  be  viewed  as  a  direct  function  of  BT.   That 
is,  for  a  constant  value  of  signal  to  noise  ratio,  the  proba- 
bility of  error  is  plotted  versus  BT.   This  plot  consists  of 
a  family  of  curves  corresponding  to  values  of  jammer  to  signal 
ratio  (J/S)  of  0.1  (-10  dB) ,  0.25  (-6  dB) .  0.5  (-3  dB)  and 
1  (0  dB) .   Finally  the  third  type  of  plot  used  in  this  analy- 
sis is  again  a  probability  of  error  versus  signal  to  noise 
ratio  plot.   However  this  plot  will  consist  of  a  family  of 
curves  for  values  of  BT  appropriate  to  the  case  analyzed.   In 
these  plots  the  value  of  J/S  used  is  zero.   This  corresponds 
to  the  case  where  the  jammer  waveform  is  considered  not  to  be 
time  truncated. 

The  case  J/S  =  0  represents  an  alternate  to  the  assumption 
that  n.(t)  is  a  time  truncated  process.   When  n.(t)  is  a  pure 
cosine  function,  N.(oo)  is  a  delta  function,  which  can  be 
completely  removed  by  the  front-end  filter.   The  resulting 
plots  then  describe  the  effect  of  the  signal  being  distorted 
by  the  front-end  filter  before  it  is  correlated  with  s,(t) . 
The  probability  of  error  curves  for  this  analysis  is  shown  by 
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the  set  of  plots  described  in  the  preceding  paragraph.  This 
performance  analysis  can  also  be  viewed  in  each  of  the  other 
plots  with  the  curve  J/S  =  0. 

B.   PHASE  SHIFT  KEY  RECEIVER  ANALYSIS 
1 .   General 

The  probability  of  error  performance  of  the  BPSK  receiver 
with  both  the  ideal  and  second  order  front-end  filters  is 
analyzed  and  compared  to  that  of  the  unfiltered  BPSK  receiver 
in  this  section.   The  results  of  Equations  3.32  and  3.33  are 
substituted  into  Equation  3.29,  yielding  the  average  proba- 
bility of  error  for  the  BPSK  receiver  with  an  ideal  front-end 
filter.   Substituting  Equations  3.69  and  3.70  into  Equation 
3.66  yields  the  average  probability  of  error  for  the  BPSK 
receiver  with  a  second  order  front-end  filter.   In  each  case 
the  plots  described  in  the  beginning  section  of  this  chapter 
are  used  to  analyze  the  performance.   For  PSK  modulated  sig- 
nals, the  discrete  values  of  BT  are  chosen  with  respect  to 
the  PSK  modulation  frequency  spectrum.    Recalling  that  the 
spectrum  (Equation  3.22)  is  of  Sa(ooT/2)  form,  the  first  null 
of  this  spectrum  occurs  for  go  =  2tt/T.   Thus  the  null  to  null 
bandwidth  is  4tt/T.   Therefore,  the  value  of  BT  =  4tt  will  be 
used  as  a  practical  limit  in  this  analysis.   The  other  dis- 
crete values  of  BT  are:   BT  =  0  (no  filter) ,  BT  =  0 . 4  (a 
small  notch),  BT  =  3.16  (25%  of  the  null  to  null  bandwidth), 
BT  =  6.28  (50%  of  the  null  to  null  bandwidth),  and  BT  =  12.56 
(the  full  null  to  null  bandwidth) . 
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2.   PSK  Receiver  With  An  Ideal  Front-End  Filter 

Figures  A. 3  through  A. 9  result  when  the  probability  of 
error  for  the  BPSK  receiver  with  a  front-end  filter  is  calcu- 
lated using  the  ideal  filter  model.   These  plots  show  how  the 
probability  of  error  performance  is  affected  by  the  signal  to 
noise  ratio,  jammer  to  signal  ratio  and  BT.   Figure  A. 3  repre- 
sents the  probability  of  error  performance  of  the  BPSK  receiver 
with  no  front-end  filter.   This  corresponds  to  the  results 
derived  and  plotted  in  Ref.  2.   In  fact  the  J/S  =  0  curve  in 
this  plot  is  the  well-known  result  for  the  coherent  BPSK  re- 
ceiver operating  in  additive  white  Gaussian  noise  alone.   This 
curve  will  serve  as  a  performance  reference  in  the  following 
discussions.   Figures  A. 3  through  A. 8  allow  the  probability  of 
error  performance  to  be  graphically  analyzed  for  the  cases 
where  the  jammer  is  either  a  truncated  process  or  a  continuous 
time  function  as  discussed  in  the  first  section  of  this  chapter. 

For  the  truncated  jammer,  by  following  the  J/S  =0.1 
and  0.5  curves  from  Figure  A. 3  to  Figure  A. 8,  the  curves  shift 
to  the  right  as  the  value  of  BT  is  increased.   This  means  that 
for  a  given  probability  of  error,  the  receiver  operating  with 
the  larger  value  of  BT  requires  more  signal  to  noise  ratio 
to  achieve  a  given  level  of  performance.   Thus  for  the  optimum 
jammer  (truncated  by  the  receiver  process)  the  front-end  filter 
actually  further  degrades  the  probability  of  error  performance 
for  the  BPSK  receiver.   Specifically,  (refer  to  the  J/S  =  0 
curve  of  Figure  A. 3)  for  BT  =  3.16  there  is  5  dB  of  degradation, 
for  BT  =  6.2  8  there  is  6  dB  of  degradation  and  for  BT  =  12.56 
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there  is  9  dB  of  degradation.   The  rapid  increase  of  the 
probability  of  error  as  BT  increases  is  seen  in  Figure  A. 8. 
For  a  fixed  value  of  E/N  =  100  (20  dB) ,  the  probability  of 
error  rapidly  increases  as  BT  approaches  4tt.   As  BT  increases 
further,  the  slope  of  the  probability  of  error  curves  become 
less  positive  as  less  energy  is  incrementally  removed  from  the 
signal  by  the  filter. 

Assuming  now  that  the  jammer  is  not  truncated  by  the 
receiver  process,  the  analysis  of  the  probability  of  error 
curves  is  performed  by  comparing  the  J/S  =  0  curves  of  Figures 
A. 4  to  A. 8  with  that  of  Figure  A. 3.   These  are  plotted  together 
in  Figure  A. 9.   From  Figure  A. 3,  for  J/S  equal  to  -10  dB , 
5  dB  of  additional  signal  to  noise  ratio  is  required  to  achieve 
the  same  probability  of  error  performance  as  with  no  jammer 
present.   For  J/S  equal  to  -3  dB ,  10  dB  of  additional  signal 
to  noise  ratio  is  required,  and  for  J/S  greater  than  one,  no 
increase  in  signal  to  noise  ratio  results  in  improved  perfor- 
mance as  discussed  in  Ref.  2.   Assuming  that  the  filter  com- 
pletely removes  the  jammer,  Figure  A. 9  shows  the  performance 
penalty  associated  with  the  use  of  the  front-end  filter.   By 
comparing  Figure  A. 9  with  Figure  A. 3,  it  is  seen  that  for 
large  values  of  J/S  the  receiver  probability  of  error  perfor- 
mance is  improved  by  the  front-end  filter.   The  specific  level 
of  improvement  is  determined  by  comparing  the  curve  in  Figure 
A. 3  corresponding  to  a  particular  J/S  value  with  the  curve  in 
Figure  A. 9  for  a  fixed  value  of  BT .   For  the  ideal  jammer  in 
which  the  frequency  spectrum  is  a  delta  function,  an  infinitely 
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narrow  bandwidth  filter  can  be  used  to  counter  this  jammer. 
In  a  more  practical  sense,  the  filter  bandwidth  should  be  made 
as  narrow  as  possible,  and  yet  completely  null  the  jammer  over 
its  entire  spectrum.   The  sensitivity  of  the  probability  of 
error  performance  to  the  value  of  BT  used  for  this  case  is 
shown  in  Figure  A. 8  with  the  J/S  =  0  curve.   Again,  there  is 
a  rapid  increase  of  probability  of  error  as  BT  increases  to 
the  value  of  12.56.   This  further  shows  that  BT  should  be  kept 
as  small  as  practically  possible. 

3.   PSK  With  A  Second  Order  Front-End  Filter 

Figures  A. 10  through  A. 16  result  when  the  probability 
of  error  for  a  coherent  BPSK  receiver  with  a  second  order 
front-end  filter  is  plotted  versus  E/N,  J/S,  and  BT.   The 
probability  of  error  in  this  case  is  obtained  by  substituting 
Equations  3.69  and  3.70  into  Equation  3.66.   Figure  A. 10 
shows  the  probability  of  error  performance  with  BT  =  0  (actually 
BT  =  0.0001  for  computational  purposes) .   For  values  of  BT 
less  than  3.16  (approximately  equal  to  tt)  ,  the  probability  of 
error  curves  compare  closely  with  the  curves  plotted  for  the 
BPSK  receiver  with  the  ideal  filter.   This  is  illustrated  by 
comparing  Figure  A. 11  with  Figure  A. 4,  and  Figure  A. 12  with 
Figure  A.  5.   For  values  of  BT  greater  than  2tt  the  comparison 
of  Figures  A. 6  and  A. 7,  with  Figures  A. 13  and  A. 14,  show  that 
for  BT  =  6.28  the  ideal  filter  calculations  differ  from  the 
second  order  filter  calculations  by  1  dB .   For  BT  =  12.56 
there  is  a  2  dB  difference  in  the  curves  calculated  using  the 
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ideal  versus  the  second  order  filter.   The  comparison  of 
Figures  A. 8  and  A. 15  show  that  as  BT  increases,  the  calculations 
using  the  ideal  filter  are  less  accurate  compared  to  the  second 
order  filter  calculations.   Again,  for  the  truncated  jammer 
case,  no  value  of  BT  results  in  improved  probability  of  error 
performance  in  comparison  to  the  unfiltered  receiver. 

Now,  assuming  that  the  receiver  process  does  not  trun- 
cate the  jammer,  Figure  A. 16  shows  the  probability  of  error 
performance  achieved  for  BT  =  0,  3.16,  6.28  and  12.56.   Com- 
paring Figure  A. 16  to  Figure  A. 9  shows  that  there  is  similarity 
in  the  probability  of  error  performance  calculated  using  the 
ideal  filter  model  with  that  calculated  using  the  second  order 
filter  model,  for  values  of  BT  less  than  3.16.   Therefore  for 
values  of  BT  less  than  it  (or  B  less  than  25%  of  the  signal  null 
to  null  bandwidth) ,  the  probability  of  error  can  be  determined 
from  the  results  derived  for  the  ideal  filter  model  being  used 
as  a  front-end  filter. 

C.   COHERENT  FREQUENCY  SHIFT  KEY  RECEIVER  ANALYSIS 

The  probability  of  error  performance  of  the  FSK  receiver 
with  the  two  filters,  the  ideal  and  second  order  front-end 
filters,  is  analyzed  and  compared  to  the  unfiltered  coherent 
FSK  receiver  performance  in  this  section.   The  results  of 
Equations  4.53  and  4.54  are  substituted  into  Equation  4.48, 
yielding  the  average  probability  of  error  of  the  FSK  receiver 
with  an  ideal  front-end  filter.   Substituting  Equations  2.28 
and  2.29  into  Equation  2.23  yields  the  average  probability  of 
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error  for  the  FSK  receiver  with  a  second  order  front-end  filter. 
In  each  case  the  plots  described  in  the  beginning  section  of 
this  chapter  are  used  to  analyze  performance.   For  FSK  modu- 
lated signals,  not  only  must  the  parameter  BT  be  chosen,  but 
also  the  frequency  separation  of  the  FSK  signals  must  be  set. 
From  the  general  results  on  performance  for  the  coherent  re- 
ceiver, the  probability  of  error  is  given  by  (from  Ref.  1 
assuming  no  filters  and  no  jamming  present) 


P    =   erfc[ — £(1  p) (5.1) 
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/NQ£(l-p) 
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where  e  and  p  have  been  defined  by  Equations  4.24  and  4.26, 
respectively.   The  results  of  Ref.  2  given  by  Equation  4.5, 
with  the  Filter  Factor  equal  to  2  ,  use  a  value  of  p  of  zero. 
Therefore  the  comparisons  made  here  will  be  for  values  of 
u),-con  such  that  p  is  zero.   The  term  used  to  describe  this 
separation  and  used  in  the  derivation  of  the  probability  of 
error  results  is  co-,T,  defined  in  Equation  4.7  as  w,T  =  (oo, -co-)  T/2 
The  minimum  value  of  oj-,T  is  chosen  to  be  tt/2  because  from 
Equation  4.27  it  is  seen  that  the  separation  (co,-con  =  tt/T)  is 
the  minimum  separation  required  for  the  normalized  correlation 
coefficient  of  sQ(t)  and  s, (t)  to  be  zero.   The  other  value  of 
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w,T  chosen  is  it,  where  again  with  this  separation  the  normalized 
correlation  coefficient  is  zero. 

In  Section  IV. A. 1  the  strategy  of  placing  the  filter  at 
the  midpoint  of  u)]_  and  oon  is  described.   From  Equation  4.13, 
it  is  seen  that  the  optimum  jammer  for  the  coherent  FSK  receiver 
consists  of  Sa(  T/2)  functions  centered  at  dii  and  w0 .   There- 
fore in  order  to  filter  out  the  major  portions  of  the  jammer, 
the  bandwidth  of  the  filter  must  be  at  least  equal  to  the  fre- 
quency separation  of  the  FSK  signals.   This  also  holds  true 
for  the  case  where  the  jammer  is  considered  not  to  be  truncated 
by  the  receiver  process.   Given  these  conditions,  for  co,T  equal 
to  tt/2  the  values  of  BT  used  will  be  zero,  corresponding  to 
no  filter,  it  corresponding  to  the  difference  of  the  signal 
frequencies,  and  finally  tt  +  O.Itt  (=  3.46)  in  order  to  take 
into  account  the  non-truncated  jammer  having  non-zero  spectral 
width.   Similarly  for  oo-jT  =  tt  ,  the  values  of  BT  chosen  are 
BT  =  0,  BT  =  2  ,  and  BT  =  2tt  +  0.1'2tt  (=  6.91). 

1.   FSK  Receiver  With  An  Ideal  Front-End  Filter 

Figures  A. 17  to  A. 24  result  when  the  probability  of 
error  for  the  coherent  FSK  receiver  with  a  single  ideal  front- 
end  filter  is  plotted  as  functions  of  signal  to  noise  ratio, 
jammer  to  signal  ratio,  frequency  separation,  and  BT .   Figures 
A. 17  and  A. 20  present  the  probability  of  error  for  the  coherent 
FSK  receiver  operating  with  no  front-end  filter.   Also  the 
J/S  =  0  curve  in  these  figures  correspond  to  the  probability 
of  error  performance  for  the  FSK  coherent  receiver  (with 
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orthogonal  signals)  operating  in  white  Guassian  noise  alone. 
It  is  not  surprising  that  Figures  A. 17  and  A. 20  are  exact 
duplicates,  because  in  each  case  the  normalized  correlation 
coefficient  of  s, (t)  and  sn(t)  is  zero,  as  described  in 
Section  V.C.   As  with  the  BPSK  analysis,  the  J/S  =  0  curves 
of  Figures  A. 17  and  A. 20  will  serve  as  a  performance  reference 

Assuming  the  truncated  version  of  the  jammer,  first 
for  oo, T  =  it/2,  Figures  A.  17  to  A.  19  show  that  as  BT  increases 
from  0  to  3.46,  the  probability  of  error  curves  shift  to  the 
right.   This  indicates  that  for  increasing  filter  bandwidths, 
more  signal  to  noise  ratio  is  required  to  achieve  the  same 
probability  of  error  performance  as  with  no  filter.   From 
Figure  A.  17,  for  a  value  of  BT  equal  to  it,  approximately  2 
dB  of  additional  signal  to  noise  ratio  is  required  to  achieve 
the  same  level  of  performance  as  in  Figure  A. 17.   This  applies 
for  all  values  of  J/S.   From  Figure  A. 19,  2.8  dB  of  additional 
signal  to  noise  ratio  is  required  to  achieve  the  performance 
level  presented  in  Figure  A. 17.   The  increase  in  probability 
of  error  with  increasing  BT  for  all  values  of  J/S  is  shown  in 
Figure  A. 23. 

With  (jo,T  equal  to  it,  more  filter  bandwidth  is  required 
in  order  to  attempt  to  notch  out  the  jammer.  As  a  result,  a 
greater  amount  of  the  signal  energy  is  removed  by  the  filter 
than  is  the  case  for  the  receiver  operating  with  oj-,T  of  tt/2. 
Therefore,  the  performance  of  the  FSK  receiver  with  a  front- 
end  filter  and  co  ,T  set  to  the  value  of  tt  ,  will  be  worse. 

d 

Figures  A. 20  to  A. 22  show  that  for  BT  increasing  from  zero  to 
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6.91,  more  signal  to  noise  ratio  is  required  to  achieve  the 

performance  level  of  the  receiver  operating  with  an  unfiltered 

jammer.   Specifically,  5.6  dB  of  additional  E/N  is  required 

for  BT  =  2tt,  and  6.6  dB  of  additional  E/N  is  required  for 

BT  =  6.91.   Figure  A. 24  shows  the  sensitivity  of  the  probability 

of  error  performance  for  increasing  values  of  BT.   Focusing 

on  Figures  A. 17  to  A. 24,  in  either  case  (co^T  =  tt/2  or  1.57), 

analysis  using  the  truncated  jammer  assumption  shows  that  no 

value  of  BT  causes  a  decrease  in  the  probability  of  error  for 

the  coherent  FSK  receiver.   Only  degradation  is  observed. 

The  analysis  of  the  coherent  FSK  receiver  with  the 

ideal  front-end  filter,  assuming  the  non-truncated  jammer,  is 

performed  by  comparing  the  J/S  =  0  curves  on  Figures  A. 17, 

A.  18  and  A.  19  for  co,T  =  tt/2.   For  w  .T  =  tt  ,  the  J/S  =  0  curves 

d  d 

on  Figures  A. 20,  A. 20  and  A. 22  are  compared.   As  before,  the 

unfiltered  J/S  =  0  curves  (Figures  A. 17  and  A. 20)  serve  as  a 

reference.   As  in  the  BPSK  case,  this  analysis  shows  that  for 

either  aj,T  =  tt/2  or  tt  ,  the  J/S  =  0  curves  still  shift  to  the 

right  as  BT  is  increased  indicating  that  more  signal  to  noise 

ratio  is  required  in  order  to  achieve  the  same  performance  as 

for  the  unfiltered  case.   From  the  curves  for  an  unfiltered 

jammer  (Figures  A. 17  and  A. 20)  for  values  of  J/S  greater  than 

0.5,  the  receiver  is  essentially  inoperable  (P  z    0.5).   However, 

the  J/S  =  0  curve  of  Figure  A. 18  shows  that  an  increase  of 

E/Nn of  5.6  dB  can  bring  the  performance  of  the  filtered  FSK 

receiver  to  the  same  value  of  P   as  for  the  unfiltered  receiver. 

e 

Figure  A.  18  corresponds  to  co,T  =  tt/2,  BT  =  tt  .   Figure  A.  19 
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shows  that  for  BT  =  3.4  6,  2.8  dB  of  additional  E/N-  is  required 
to  achieve  the  unfiltered,  un jammed  performance.   Figures  A. 18 
and  A. 19  show  that  for  values  of  J/S  greater  than  .5  (for  the 
coherent  receiver)  the  use  of  the  ideal  filter  does  show  an 
improvement  on  the  performance  of  the  FSK  receiver  when  the 
filter  is  used  to  notch  out  the  cw  jammer.   The  sensitivity 
of  the  probability  of  error  to  the  value  of  BT  used  for  the 
case  of  oo, T  =  tt/2,  is  shown  by  the  J/S  =  0  curve  of  Figure 
A. 23.   As  in  the  BPSK  receiver  case,  this  shows  that  the  smallest 
possible  filter  bandwidth  should  be  utilized.   For  (jo,T  =  tt  , 
similar  results  are  shown  in  Figures  A. 21  and  A. 22.   With 
BT  =  2tt,  5.6  dB  of  additional  E/N_  is  required  to  achieve  the 
performance  level  of  the  unfiltered,  un jammed  receiver,  and 
for  BT  =  6.91,  6.6  dB  of  additional  E/Nn  is  required.   Figure 
A. 24  shows  the  sensitivity  of  the  probability  of  error  for 
this  case.   An  important  fact  to  remember  is  that  when  using 
Figures  A. 23  and  A. 24,  the  ideal  filter  does  not  alter  the 

jammer  until  BT  is  at  least  2xco,T. 

J  d 

2 .   The  FSK  Coherent  Receiver  With  A  Second  Order 
Front-End  Filter 

Figures  A. 25  to  A. 32  result  when  Equations  4.53  and 

4.54  are  substituted  into  Equation  4.48  and  are  plotted  versus 

BT,  J/S  and  E/N  .   Figures  A. 25  and  A. 28  are  the  probability 

of  error  curves  for  the  unfiltered  FSK  coherent  receiver 

operating  with  oj,T  =  tt/2  and  tt  ,  respectively,  and  with  an 

unfiltered  jammer.   As  described  in  the  analysis  of  Figures 

A. 17  and  A. 20,  Figures  A. 25  and  A. 28  are  exact  duplicates 
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because  p  in  each  case  is  zero.   In  general,  the  curves  plotted 
using  the  second  order  filter  model  match  within  1  dB  to  those 
calculated  using  the  ideal  filter  model.   Specifically,  Figures 
A.  26  and  A.  27  (os,T  =  tt/2,   BT  =  tt  and  3.46  respectively) 
compare  within  1  dB  to  the  results  of  Figures  A. 18  and  A. 19 
which  are  calculated  using  the  ideal  filter  model.   Similarly, 
for  w,T  =  tt  ,  Figures  A.  29  and  A.  30  compare  within  1  dB  to  the 
results  presented  in  Figures  A. 21  and  A. 22  respectively.   The 
sensitivity  analysis  from  Figures  A.  31  and  A.  32  (w,T  =  tt/2 
and  tt)  still  shows  the  rapid  increase  in  probability  of  error 
with  increasing  BT.   However,  compared  to  similar  curves  calcu- 
lated using  the  ideal  filter,  Figures  A. 31  and  A. 32  have 
smaller  (positive)  slope.   In  all,  assuming  the  truncated 
jammer,  Figures  A. 25  to  A. 32  still  show  that  as  the  front-end 
filter  bandwidth  increases,  the  probability  of  error  increases 
also.   Therefore,  as  with  the  ideal  filter  case,  the  second 
order  front-end  filter  only  further  degrades  the  performance  of 
the  coherent  FSK  receiver  operating  against  the  truncated 
optimum  jammer. 

Assuming  the  non-truncated  version  of  the  optimum 
jammer,  it  is  easily  seen  that  a  second  order  front-end  filter 
such  as  the  one  given  by  Equation  3.37  would  not  be  completely 
effective  against  that  jammer.   By  expanding  Equation  4.3, 
it  is  seen  that  the  jammer  spectrum  would  consist  of  delta 
functions  centered  at  co,  and  oon .   The  filter  described  in 

Section  III.B  has  a  null  only  at  to  =  ±co  .   Therefore  the 

J  c 
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analysis  cannot  be  performed  by  simply  letting  J/S  =  0  in 
this  case.   From  the  analysis  completed  for  the  ideal  filter 
case,  it  is  seen  that  a  higher  order  filter,  or  a  filter  with 
nulls  at  00-,  and  con  ,  should  be  used  to  defeat  the  non- truncated 
optimum  jammer. 
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VI.   CONCLUSIONS 

A.   GENERAL  CONCLUSIONS 

In  this  thesis  the  effectiveness  of  front-end  filtering 
techniques  hasi  been  investigated  as  an  electronic  counter- 
countermeasure  (ECCM)  in  order  for  the  coherent  digital 
receiver  to  be  able  to  operate  effectively  against  the  optimum 
jammer  derived  in  Ref.  2.   The  specific  coherent  digital  re- 
ceivers analyzed  were  the  coherent  BPSK  receiver  and  the 
coherent  BFSK  receiver.   The  filters  placed  in  the  front-end 
of  the  receivers,  as  shown  in  Figure  A. 2,  were  an  ideal  filter 
and  a  second  order  filter  of  the  form  given  by  Equation  3.37. 

Because  the  analysis  of  Ref.  2  defines  the  optimum  jammer 
in  the  time  interval  (0,T) ,  the  analysis  was  performed  using 
this  interval  in  the  calculation  of  the  Fourier  transforms 
of  the  jammer  models.   As  a  result,  Sections  V.B.2  and  V.B.3 
as  well  as  Sections  V.C.2  and  V.C.3  show  that  none  of  the 
front-end  filters  resulted  in  any  decrease  in  the  probability 
of  error  in  comparison  to  receivers  operating  without  front- 
end  filters  for  the  class  of  receivers  analyzed.   This  occurs 
because  the  resultant  jammer  consists  of  a  signal  that  is  of 
the  same  form  as  s,  (t)  ,  the  correlation  signal.   As  a  result, 
the  filter  affects  equal  portions  of  both  the  signal  spectrum 
as  well  as  the  jammer  spectrum.   Therefore  under  this  analysis, 
the  front-end  filters  further  added  to  the  probability  of  error 
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caused  by  the  jammer  and  the  additive  noise  for  both  the  BPSK 

and  BFSK  coherent  receivers. 

Taking  into  account  the  physical  aspect  of  this  problem, 

an  ambiguity  is  seen.   Both  s,(t)  and  n.(t)  are  defined  over 

d         j 

the  interval  (0,T).   However  over  the  next  interval  (T,2T) 

these  signals  remain  unchanged.   This  repeats  itself  for  the 

next  interval  of  length  T,  and  so  on.   In  other  words,  both 

s-,(t)  and  n.(t)  are  cw  signals.   Assuming  such  a  form  for 
J 

the  optimum  jammer,  the  analysis  was  carried  out  by  substi- 
tuting 0  for  the  value  of  J/S ,  for  all  but  the  FSK  receiver 
with  the  second  order  front-end  filter.   Setting  J/S  =  0  pro- 
vides the  same  result  as  setting  the  term  (n!,s,)  equal  to 
zero.   In  the  case  of  the  FSK  coherent  receiver,  the  second 
order  filter  does  not  have  nulls  at  the  frequencies  of  the 
optimum  jammer.   This  analysis  cannot  be  performed  simply  by 
setting  J/S  =  0 . 

When  the  PSK  coherent  receiver  was  analyzed  with  both 
filter  models  using  the  cw  form  of  the  optimum  jammer,  the 
filter  was  observed  to  introduce  improvements  to  the  proba- 
bility of  error  performance  provided  that  the  receiver  was 
jammed  with  a  value  of  J/S  of  at  least  0.5.   This  is  quite 
apparent  for  J/S  equal  to  1.0,  as  the  unfiltered  PSK  coherent 
receiver  becomes  essentially  inoperative  (P  =  0.5).   When 
the  filter  is  introduced  in  this  case,  the  jammer  term 

(n!,s  )  is  set  to  zero.   Therefore  the  conditional  mean  of 
J   d 

G*,  the  output  of  the  receiver,  is  a  function  of  (s!,s,)  only. 
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When  this  happens,  any  increase  in  the  probability  of  error 
compared  to  that  of  the  unfiltered,  un jammed  receiver  is 
caused  only  by  the  filter  distorting  the  signal.   The  use  of 
the  filter  is  not  without  cost.   As  shown  in  Figures  A. 9  and 
A. 16,  additional  signal  to  noise  ratio  is  required  to  achieve 
the  same  probability  of  error  as  that  of  unfiltered  unjammed 
receiver. 

The  analysis  of  the  FSK  coherent  receiver  with  the  ideal 
front-end  filter,  under  the  assumption  of  the  cw  jammer  shows 
similar  results.   However,  because  a  single  filter,  centered 
at  the  midpoint  frequency  between  the  FSK  signals  is  used, 
more  signal  to  noise  ratio  is  required  to  achieve  the  same 
probability  of  error  performance  as  that  of  unfiltered,  un- 
jammed receiver.   This  results  because  the  filter  bandwidth 
must  be  at  least  equal  to  the  frequency  separation  of  the  FSK 
signals.   Due  to  this  observation,  the  smallest  practical  fre- 
quency separation  should  be  used.   Analysis  of  the  second  order 
front-end  filter  shows  that  this  type  of  filter  centered  midway 
between  the  FSK  signals  was  not  as  useful  as  a  filter  with 
nulls  at  the  jammer  frequencies. 

B.   SPECIFIC  CONCLUSIONS 

Some  specific  conclusions  follow  from  the  analysis  per- 
formed in  this  thesis.   These  are  as  follows: 

1)   For  the  PSK  coherent  receiver  with  a  front-end  filter, 
the  probability  of  error  performance  calculated  using 
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the  ideal  front-end  filter  matches  the  performance 
calculated  using  the  second  order  filter  within  1  dB 
for  values  of  BT  of  less  than  tt.   Therefore,  the  analysis 
can  be  performed  using  the  ideal  filter  for  small  values 
of  BT.   So,  though  this  will  require  computer  integra- 
tion methods,  time  will  be  saved  by  not  analyzing  complex 
filter  models. 

2)  Under  the  assumption  of  a  time  truncated  jammer,  the 
probability  of  error  for  the  FSK  coherent  receiver 
calculated  using  the  ideal  filter  compares  with  the 
results  calculated  using  the  second  order  filter  with 
BT  =  2to,T.   Therefore,  the  ideal  filter  can  be  used  for 
these  calculations. 

3)  The  second  order  front-end  filter  centered  between 
the  FSK  signals  does  not  effectively  remove  the  jammer 
(assuming  the  non-truncated  jammer)  and  therefore  is 
not  as  useful  as  a  filter  with  nulls  at  the  jammer 
frequencies . 

4)  It  appears  that  part  of  the  performance  degradation 
introduced  by  filtering  is  due  to  the  fact  that  the 
filter  distorts  the  signal  components  s, (t)  and  sn(t), 
so  that  the  correlator  receiver  (or  equivalently  the 
matched  filter  receiver)  is  not  matched  to  the  incoming 
signals.   This  degradation  could  be  overcome  by  using  a 
matched  filter  receiver  that  is  matched  to  s* (t)  and 

s' (t) .   While  this  has  not  been  investigated  in  this 
thesis,  it  appears  that  such  investigations  would  be  of 
great  interest  and  worthwhile. 
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Figure   A.  3.       Coherent   PSK  Receiver  with    an    Ideal 
Front-End    Filter,    BT    =    0.0 
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Figure   A. 4 


Coherent   PSK  Receiver  with   an    Ideal 
Front-End    Filter,    BT   =    0.4 
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Figure  A. 5.   Coherent  PSK  Receiver  with  an  Ideal 
Front-End  Filter,  BT  =  3.16 
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Figure  A. 6.   Coherent  PSK  Receiver  with  an  Ideal 
Front-End  Filter,  BT  =  6.28 
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Figure  A. 7.   Coherent  PSK  Receiver  with  an  Ideal 
Front-End  Filter,  BT  =  12.56 
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Figure   A. 8. 


Coherent  PSK  Receiver  with  an  Ideal 
Front-End  Filter,  E/N  =  100 
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Figure   A. 9 


Coherent    PSK    Receiver   with    an    Ideal 
Front-End   Filter  with   Variable   BT 
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Figure  A. 10.   Coherent  PSK  Reciver  with  a  Second 
Order  Front-End  Filter  (w/BT  =  0.0) 
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Figure  A. 11 


Coherent  PSK  Receiver  with  a  Second 
Order  Front-End  Filter  (w/BT  =0.4) 
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Figure   A. 12. 


Coherent  PSK  Reciver  with  a  Second 
Order  Front-End  Filter  (w/BT  =  3.16) 
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Figure   A. 13. 


Coherent    PSK   Receiver  with    a   Second 
Order   Front-End   Filter    (w/BT   =    6.28) 
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Figure   A. 14. 


Coherent  PSK  Receiver  with  a  Second 
Order  Front-End  Filter  (w/BT  =  12.56) 
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Figure  A. 15 


Coherent  PSK  Receiver  with  a  Second 
Order  Front-End  Filter  (E/N  =  100) 
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Figure  A. 16.   Coherent  PSK  Receiver  with  a  Second 

Order  Front-End  Filter  (w/various  BT) 
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Figure  A. 17.   Coherent  FSK  Receiver  with  an  Ideal 

Front-End  Filter,  WdT  =  1.57,  BT  =  0.0  0 
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Figure  A. 18.   Coherent  FSK  Receiver  with  an  Ideal 

Front-End  Filter,  WdT  =  1.57,  BT  =  3.14 
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Figure  A. 19.   Coherent  FSK  Receiver  with  an  Ideal 

Front-End  Filter,  WdT  =  1.57,  BT  =  3.4  6 
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Figure   A. 20 


Coherent  FSK  Receiver  with  an  Ideal 
Front-End  Filter,  WdT  =  3.14,  BT  =  0.00 
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Coherent    FSK   Receiver  with   an    Ideal 
Front-End   Filter,    WdT   =    3.14,    BT   =    6.2  8 
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Figure  A. 22.   Coherent  FSK  Receiver  with  an  Ideal 

Front-End  Filter,  WdT  =  3.14,  BT  =  6.91 
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Figure   A. 23. 


Coherent   FSK   Receiver  with   an    Ideal 
Front-End   Filter,    WdT   =    1.57,    E/N   =10  0 
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Figure   A. 24. 


Coherent  FSK  Receiver  with  an  Ideal 
Front-End  Filter,  WdT  =  3.14,  E/N  =  100 
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Figure   A. 25.       Coherent   FSK   Receiver  with   a   Second 

Order   Front-End    Filter,    WdT   =    1.57,    BT   =    0.00 
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Figure  A. 26.   Coherent  FSK  Receiver  with  Second  Order 
Front-End  Filter,  WdT  =  1.57,  BT  =  3.14 
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Figure  A. 27. 


Coherent  FSK  Receiver  with  a  Second 

Order  Front-End  Filter,  WdT  =  1.57,  BT  =  3.46 
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Coherent  FSK  Receiver  with  a  Second  Order 
Front-End  Filter,  WdT  =  3.14,  BT  =  0.00 
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Figure  A. 29 


Coherent  FSK  Receiver  with  a  Second  Order 
Front-End  Filter,  WdT  =  3.14,  BT  =  6.28 
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Coherent  FSK  Receiver  with  a  Second  Order 
Front-End  Filter,  WdT  =  3.14,  BT  =  6.91 
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Coherent  FSK  Receiver  with  a  Second  Order 
Front-End  Filter,  WdT  =  1.57,  E/N  =  100 
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Figure   A. 32 


Coherent  FSK  Receiver  with  a  Second  Order 
Front-End  Filter,  WdT  =  3.14,  E/N  =  100 
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